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ABSTRACT 


The optical path length from a satellite to the earth's surface Is 
strongly dependent on the atmospheric pressure along the propagation path. 
Surface pressure can be determined by measuring the difference between the 
round-trip propagation times of laser pulses that are transmitted simulta- 
neously at two wavelengths. Although pressure measurements can be made 
over the ground and water, the application of this technique to pressure 
measurements over the ocean Is considered. 

The statistical characteristics and the waveforms of the ocean- 
reflected laser pulses are studied. The received signal Is found to be 
corrupted by shot noise and time-resolved speckle. The statistics of time- 
resolved speckle and Its effects on the timing accuracy of the receiver are 
studied In the general context of laser altimetry. 

For estimating the differential propagation time, various receiver 
timing algorithms are proposed and their performances evaluated. The 
results Indicate that, with the parameters of a realistic altimeter, a 
pressure measurement accuracy of a few millibars Is feasible. 

The data obtained from the first airborne two-color laser altimeter 
experiment are processed and analyzed. The results are used to verify the 
pressure measurement concept. 
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1 . INTRODUCTION 

Global measurements of temperature and pressure are essential In all 
weather and climate predictions. Forecasting and modeling have been ham- 
pered by the lack of measurements In Inaccessible regions of the earth's 
surface, particularly over the oceans. Consequently, there Is considerable 
Interest In developing remote-sensing techniques for measuring pressure and 
temperature from satellites. Pressure Is probably the most difficult para- 
meter to measure accurately. Korb et al. [1] have proposed a differential 
absorption (DIAL) technique for remote sensing of pressure and temperature 
by observing the near-IR absorption in the oxygen A band. The feasibility 
of their technique was demonstrated experimentally using CW lasers to make 
measurements over the horizontal paths. Unfortunately, the technique 
places severe constraints on the wavelength stability of the laser. 

In a previous paper, Gardner [21 proposed a pressure measurement tech- 
nique using a two-color short pulse laser altimeter. This technique was an' 
extension of NASA's ongoing two-color ranging program to measu^re atmo- 
spheric pressure [3]. The technique makes use of the fact that the 
atmosphere Is dispersive and the group refractlvlty Is proportional to 
pressure; therefore, the difference between the optical path lengths from a 
satellite to the earth's surface Is proportional to the surface pressure. 
The altimeter uses fixed frequencies rather than tunable lasers and Is 
Insensitive to laser wavelength changes. The theory on which the technique 
Is based Is reviewed In Chapter 2. 

Although pressure measurements can be made over both the ground and 
water, pressure measurement over the ocean Is of greatest Interest. In a 
previous paper [4], we studied the statistics and waveforms of the ocean 
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reflected laser pulses. Parts of the results that are relevant to the pres- 
sure measurement technique are reviewed In Chapter 3. In addition, we also 
study the reflection of laser pulses from sinusoidal and trocholdal waves, 
because these waves give rise to glints that can Improve the accuracies of 
the pressure measurements. 

When a short laser pulse Is reflected from the ocean surface, the 
reflected pulse will be broadened to about twice the width of the range 
spread of the wave height. Because the received pulse has a width longer 
than the transmitted laser pulse, speckle will cause random small-scale 
fluctuations within the received pulse which distort Its shape. This 
effect Is referred to as time-resolved speckle. In Chapter 4, we derive 
the statistics of time-resolved speckle. The problem of estimating the 
arrival time of laser pulses In the presence of time-resolved speckle Is 
also considered. The results not only are Important to the pressure 
measurement technique, but also have applications In general laser radar 
and ranging. 

The pressure measurement technique requires accurate measurements of 
the differential propagation time of laser pulses at two wavelengths. If 
the altimeter uses the fundamental and tripled ND:YAG laser frequencies, 
picosecond timing accuracies are required to make mllllbar-level pressure 
measurements. The timing problem Is complicated by the fact that, due to 
the random dynamic nature of the ocean, the mean reflected pulse shape can- 
not be predicted a priori. This means optimal estimators such as the maxi- 
mum likelihood estimator cannot be used. In Chapter 5, we propose and 
study several suboptlmal estimators that do not require knowledge of the 
pulse shape. 

The timing accuracy of the receiver depends on the bandwidth of the 
received signal. The signal bandwidth Is related to the surface profile 
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within the laser footprint. In Chapter 6, we investigate the relation be- 
tween the surface profile and the bandwidth of the received signal. 

In Chapter 7, we use the results obtained from previous chapters to 
estimate the expected timing accuracies of realistic system designs and to 
estimate the corresponding accuracies of the barometric measurements. 

An airborne altimeter experiment was conducted by NASA personnel on 
September 7 and 8, 1983. This initial experiment was basically a shakedown 
flight for a series of more extensive flights to be conducted during the 
late summer 1984. To verify the theory, the data collected were processed 
by different timing algorithms. The results are discussed in Chapter 8. 
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2. THEORY OF THE PRESSURE MEASUREMENT TECHNIQUE 


The optical path length is defined as the integral of the group 
refractive index along the ray path. Because the horizontal refractivity 
gradients in the atmosphere are small, the two-way optical path length be- 
tween a satellite and the earth's surface for a pulsed laser system is 
given by [2] 


r_at a + 

R - 2 / dr 

o ^ 

r 

ocean 


sin 9 


( 2 . 1 ) 


N is the group refractivity, r is the geocentric altitude, and 9 is given 
S 

by Snell's law for a spherically stratified medium. It is convenient to 
express Eq. (2.1) as the sum of the straight-line path length R^ and an 
atmospheric .correction AC: 


R^ = 2Rj^ + AC 


( 2 . 2 ) 


where 


AC 


2 

r 

ocean 


sin 9 


+ 2 



ocean 


dr 

sin 9 



(2.3) 


The first term comprising AC is the velocity correction, while the second 
terra is the difference between the geometric lengths of the ray and 
straight-line paths. 

The atmospheric correction can be evaluated by using an appropriate 
model for the group refractivity. Although AC depends on the atmospheric 
pressure, temperature, and humidity along the propagation path, it Is most 
sensitive to pressure. In fact, when the laser is pointed at nadir, AC is 


5 


approximately proportional to the atmospheric pressure at the earth's sur- 
face. Therefore, surface pressure can be computed from measurements of AC. 

A pulsed laser altimeter measures the round-trip optical path length 
R^. If the straight-line distance (R^) between the satellite and the laser 
footprint on the earth's surface is known, AC can be calculated by using 
Eq. (2.2). Because R^ must be known to within a few centimeters for the 
approach to be effective, it is probably not practical. As an alternative, 
a multicolor altimeter can be used to determine AC by calculating the 
difference between the round-trip optical path lengths measured at two dif- 
ferent wavelengths: 

AR = Rj - R 2 = ACj - AC 2 . (2.4) 

A theoretical expression for AR can be obtained by evaluating the integrals 
in Eq. (2,3). This is done by using the perfect gas law, law of partial 
pressure, and the hydrostatic equation to obtain a suitable refractivity 
profile. Although the dominant variation of AR with respect to pressure is 
linear, there is a small quadratic variation which cannot be neglected. If 
Eqs. (2.3) and (2.4) are evaluated and then solved for atmospheric pressure 
at the earth's surface, we obtain [2] 

p = , (2.5) 


a = 4.73 X 10 


-8 2 

T sin^ E 3 ’ 


( 2 . 6 ) 


b = -2.357 . 10-3 


2 - K ’ 


tan^ E 


sin^ E 


(2.7) 
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d 


F( 9.H) AR sin E 
2[f(X^) - fCX^)] 


2.24 


X 10 ^ e 


(2o8) 


F(9,H) = 1 + 0,0026 cos (20) - 0.0003H 


(2,9) 


K =■ 1,163 + 0,00968 cos (20) - 0,00104T + 0,00001435? , (2.10) 


. 0.0164 , 0.000228 
f(X) = 0.9650 + a — + r , 

r X 

where 


( 2 . 11 ) 


X = laser wavelength ( im) , 

e •> water vapor pressure at the laser footprint (mbar), 

P = atmospheric pressure at the laser footprint (mbar), 

T =• temperature at the laser footprint (®K), 

0 a colatitude of the laser footprint, 

H = altitude of the laser footprint above sea level (km), and 
E = satellite elevation angle. 


The function f(X) accounts for atmospheric dispersion which Is responsible 
for the difference In the measured pathlengths at the two wavelengths. 

To compute the surface pressure It Is necessary to know AR, E, e and 
T. Fortunately, only a crude estimate of the surface temperature Is 
required. An accuracy of 20-30® Is easy to obtain and should be adequate. 

A temperature error of 20®C would contribute less than a few tenths of a 
millibar to the pressure error for an elevation angle above 30®. Errors In 
the measured values of E are also only significant at the lower elevation 
angles. Above 50® elevation, the errors In E can be neglected provided 
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they are of the order of a millirad or less. The pressure sensitivity to 
errors in water vapor is constant with respect to elevation angle. A 
10 mbar error in water vapor pressure will contribute ~l mbar to the sur- 
face pressure error [2]. Since water vapor pressure can approach 40 mbar 
when the surface temperature and relative humidity are high, its effect 
cannot be ignored. Currently, satellite-based microwave sensors can pro- 
vide water vapor pressure information accurate to about 10% [5]. For a 
maximum water vapor pressure of 40 mbar, a 10% error in water vapor 
pressure measurement corresponds to about 0.4 mbar error in surface 
pressure error. So, by using the information collected by the microwave 
sensor or by Installing a microwave sensor together with the laser 
altimeter, the error due to water vapor can be reduced to an insignificant 
amount. 

In an actual system, the dominant error source is likely to be the 
differencial path-length measurement. The pressure sensitivity to differen- 
tial path- length errors [2] is given by 


3P ^ 0.212 sin E 
3AR "" f(X^) - fCX^) 


(mbar/mm) 


( 2 . 12 ) 


This sensitivity depends both on the choice of wavelengths and on the ele- 
vation angle. At low-elevation angles the laser pulses travel through 
more of the atmosphere so Che pressure effects are more significant. In 
addition, AR is greater for larger wavelength differences because of 
greater differences in group ref ractlvity. 3P/ 3AR is plotted vs. elevation 
angle in Fig. 2.1 for three possible wavelength combinations of fundamen- 
tal (1064 nm) , doubled (532 nm) , and tripled (355 nm) ND:YAG laser 
frequencies. When ranging at nadir using the 355 and 1064 nm wavelengths, 
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the required differential path-length accuracy is 0.6 mm/tnbar, which 
requires a timing accuracy of 2 psec/mbar. 

Although the pressure sensitivity to errors in AR decreases as the 
elevation angle is decreased, other error factors become more significant. 
The sensitivity to temperature and elevation-angle errors and the differen- 
tial path-length errors Increase as the elevation angle is decreased. 
Consequently, the best performance will be obtained with the laser pointed 
at nadir. When the laser is pointed to within a few degrees of nadir, 

Eq. (2.5) can be ..approximated by 


= -7 “ 2.12 


X 10 


2 F(6,H) sin E 
f(Xj) - f(X2> 


AR - 0.095e 


(2.13) 


In this case, the rms pressure error is approximately 



1/2 

> 


(2.14) 


where a^, o^, and are the rms errors in AR, e, and E, respectively. 
In Eq. (2.14), 


= 0.095 (mbar/mbar) , 


(2.15) 


12. s io~^P 

3E ~ tan E 


(mbar/mrad) ; 


-(2,16) 


and (3P)/(3AR) is given by Eq. (2.12) 


10 


3. LASER 8ACKSCATTER FROM THE OCEAN 

Although pressure measurements can be made over both the ground and 
water, pressure measurement over the ocean is of greatest Interest. In a 
previous paper [4], we studied the statistics and waveforms of the ocean 
reflected laser pulses. Parts of the results that are relevant to the 
pressure measurement technique are reviewed in this chapter. In addition, 
we also study the reflection of laser pulses from sinusoidal and trocholdal 
waves, because these waves give rise to glints that can Improve the accu- 
racies of the pressure measurement. 

For near normal incidence, the reflection of laser pulses from the 
ocean is due mainly to scattering by randomly distributed specular points 
on the surface. These specular points arise from ocean-wave structures 
with slopes oriented to reflect the signal back to the receiving telescope. 
Ocean waves are generated by the joint actions of wind as a disturbing 
force, and gravity and surface tension as restoring forces. The larger 
waves, with wavelengths greater than ~2 cm, are gravity waves. The smaller 
waves, for which surface tension is the restoring force, are capillary 
waves. The wavelengths of capillary waves are bounded at the small-scale 
end by viscous dissipation to mm. Since the wavelengths of even the 
smallest capillary waves are much larger than an optical wavelength, 
diffraction is not important, and the laser backscatter can be analyzed 
using geometric optics. 

If the laser beam is pointed at nadir and the reflected pulses are 
detected by a direct-detection receiver, the mean received signal con- 
ditioned on a given ocean surface profile is given by [4], [6J 
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E[S(t)lcl = S(t) = <N> / d% b2<^.z) lf(t - * h(t) , 

where 

bjj(_p,z) = |a(_p,z)|’^ S"^^(_£)// d^^|a(_p,z) 1" , 

and 

If, = + .ei - 

^ c cz c * 

Here 

<N> = expected number of detected signal photons/pulse, 

_p = (x,y) => horizontal coordinate vector on the ocean surface 
measured from the center of the footprint, 
a(_p, z) = complex amplitude cross section of the laser footprint, 

3^(_p) = power reflection coefficient of the ocean surface, 

C(_o) = ocean surface profile, 

z = altitude of the laser altimeter, 
c = velocity of light, 
f(t) = transmitted pulse amplitude, and 
h(t) = Impulse response of the receiver electronics. 

In this equation ip is the delay of the reflected pulse. The first term in 
Eq. (3.3) is the nominal distance of the target, the second term is the 
additional delay due to the curvature of the laser wavefront and the last 
terra is from the range spread of the target. The power- re fleet ion coef- 
ficient depends on the angle of incidence and sea state. For divergence 
angles of a few mllliradians or less, B^(^) is essentially constant within 


(3.1) 

(3.2) 

(3.3) 


the footprint. 
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The covariance function of the received signal, which describes the 
fluctuation of the received signal due to shot noise and speckle, is given 
by [41, [6] 

OD 

“ <N> / b2 (_P,z) / drjf(T - h(t^ - t) h(t2 - t) 

^ oo 

2 

+ ! d^_p b^(_p,z) g(tj - i{)) g(t 2 - '!») , (3.4) 

where 

Ap [ / d^_p I a(_p, z ) 1 ^ 2 (_£) ]^ 

K = -4-5 5 4-5 (3.5) 

X“^z^ / d\p|a(_p,z)| Bp(^) 

and 

g(t) = [f(t)l^ * h(t) . (3.6) 

In Eq. (3.4) K is the ratio of the receiver aperture area (A ) to the 

R 

speckle correlation area and g(t) is the point target response of the 
system. The first term on the right-hand side of Eq. (3.4) is due to shot 
noise, while the second terra is due to speckle noise. For a Gaussian laser 
cross section, K is given by [4], [6] 

/2 tan 

IT 

where 9^ is the beam divergence (half-angle measured at exp(-l/2) point). 
For typical laser altimeter configurations, K is of the order of a few 
thousand or larger. 

The wave height at any given point on the ocean is the result of many 
wave components that have been generated by the wind in different regions 
and have propagated to the point of observation. Since the motions of dif- 
ferent wave components are weakly correlated, the central limit theorem can 




(3.7) 
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be used to argue Chat Che wave height (l.e., surface profile) is Gaussian 
distributed. Therefore, over areas which are large compared Co Che longest 
waves on Che ocean, Is usually assumed to obey Gaussian statistics. 

Experimental data tend to support this hypothesis. Consequently, the wave 
slopes are also Gaussian distributed and independent of Che wave height 
[4J. Under these conditions, the power-reflection coefficient for normal 
incidence is given by [41, [7] 


6 . 


R(Q) 


^ 47t(S^ + 2 Can^ 9^) 


(3.8) 


where R(0) is the Fresnel backscatter coefficient and S is the mean-square 
(MS) value of the total slope. The factor 4 it is needed Co convert scat- 
tering cross section to reflection coefficient. An empirical relationship 
2 

between S and wind speed was derived by Cox and Munk [8] ; 

a 


a > =. 0.003 + 0.00512W , 

\3x/ \3y/ ’ 


(3.9) 


where W is the average wind speed in m/sec measured at 12.5 m above the 
mean sea level. 

The expected received pulse can be calculated by taking Che 

expectation of Eq. (3.1) over Che probability distribution of Che surface 

profile 5* In Figs. 3.1 and 3.2, we show the mean received waveforms for a 

laser altimeter that has a Gauss lan-shape point-target response g(t). The 

point-target response is assumed to have a rras range width (ca ) of 1 cm. 

S 

In Fig. 3.2, the mean waveforms are plotted for five different sea states. 
These figures show that the larger the significant wave height (SWH), the 
broader the return pulse. In Fig. 3.2, the sea state is fixed while the 
beam divergence is varied. As Che beam divergence Increases, the trailing 
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NORMALIZED 


TIME 



Figure 3.1. Mean received waveforms of the laser altimeter for different 

wind speed and Slffl. Gaussian ocean surface, normal incidence. 
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Figure 3.2. Mean received waveforms for laser altimeters with different 
beam divergence angles. Gaussian ocean surface, normal 
incidence. 
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edge of the pulse is broadened by Che wavefront curvature effects, and the 
return pulse becomes asymmetric. To obtain the sharpest return pulses, it 
is necessary to have a small beam divergence which can be achieved rela- 
tively easily with laser altimeters. This is more difficult Co achieve 
with radar altimeters due to their limitation in Che antenna size. 

The pulse shapes computed using the above approach will accurately 
predict the actual pulse shapes whenever the laser footprint is large com- 
pared to Che periods of the waves. However, if Che footprint is small, 
the local probability distribution of the surface profile within the 
footprint may be considerably different from Gaussian. For pressure- 
measurement applications, a small beam divergence angle is desired to mini- 
mize pulse broadening due to beam-curvature effects. As a consequence, 
the size of the footprint will be small and may not be significantly 
larger chan the wavelength of the long-period ocean waves. Sea-wave 
records often reveal a sinusoidal or trocholdal profile for the dominant 
long-period waves which may extend over distances of a few hundred meters. 
In these areas, the surface profile is more suitably modeled as a large 
amplitude sinusoidal (or trocholdal) wave with superimposed small amplitude 
disturbances (capillary waves or small gravity waves) of short correlation 
length. For a sinusoidal wave model, we have 

5(^) = ?g(^) + Ci(^) 

= Acos^-^ X + ej + 5j(_o) , (3.10) 

where A is the amplitude of the sinusoidal wave, A is its wavelength, 6 
is the phase angle and lj^(_p) is the small-scale disturbance. Using this 
model, Che mean received signal can be expressed as 
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E{S(t)} = / p(5^) E{S(t)l5^} 

•00 

= <»> / d% / d?j p<{,) g(t-|5.-4-b|?„d-f <i) . 

• 00 ’ 

(3.11) 


The integration over can be viewed as a convolution, so Eq. (3.11) can 
also be written as 


E{S(t)} =yE{S(t)|5i 




(3.12) 


This equation implies that we need only evaluate the mean waveform for the 
equals zero case, l.e., E {S( t) j Cj = 0}. The mean received waveform for 
any distribution of can be obalned by convolving the probability density 
function of with E{S(t)l?j^ =0}. 

Expected received pulse shapes were calculated for several values of 
beam divergence using the sinusoidal model. The results which are plotted 
in Fig. 3.3 show strong reflections from the wave crests and troughs. 

These reflections or glints occur because of the relatively large surface 
area at the crest and trough altitudes. The sharpness of the observed 
glints is directly related to the beam divergence. When the divergence is 
large, beam-curvature effects broaden the reflected pulses. The rise time 
of the crest reflection is limited only by the transmitted pulse width, 
beam divergence, and the receiver bandwidth. The tail at the trough 
reflections is primarily due to beam curvature effects. In Fig. 3.4, we 
calculated the results for three different ocean-wave heights. The wave 
height of the ocean can be Inferred from the received pulse by measuring 
the separation between the crest and trough returns. 

When the laser is pointed slightly off nadir in a direction normal to 
the wave fronts, the received pulse broadens, and multiple glints appear 
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Altitude =400 km 

Transmitted Pulse Width = I cm (33 ps) 


Wave Heights 2m 
(Trough -Crest) 



Figure 3.3. Mean received pulse shape for reflection from an ocean wave 
with a height of 2 m (trough-crest). The laser altimeter is 
at an altitude of 400 km and the beam divergence is varied 
from 100 prad to 900 prad. 
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Altitude = 400 km 

Transmitted Pulse Width = lcm(33ps) 
Beam Divergence = lOO/irod 



RANGE (m) [|(t-§^)] 

Figure 3.4. Mean received pulse shape for reflection from ocean waves with 
different wave heights. Normal incidence. 
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because of reflections from successive wave crests and troughs. This Is 
Illustrated In Fig. 3.5 for a nadir angle of 1“. A sinusoidal wave model 
was used with a wave height of 1 m (crest to trough) and a wave period of 
10 m. Crest reflections are responsible for sharp peaks In the first half 
of the pulse, while trough reflections are responsible for the structure In 
the second half. At 1® off nadir for a 10-m wave period, the range dif- 
ference between successive crests (or troughs) Is ~17.5 cm. These results 
suggest that It may be possible to measure the crest-to-t rough wave height 
and wave period as well as SWH using a short-pulse laser altimeter. This 
Is not possible using microwave altimeters because Che curvature effects 
due to the large beamwldths (~l-2®) obscure the crest and trough reflections. 

In the actual ocean, the wave crests tend to be relatively high and 
sharp, while the wave troughs are comparably smooth and shallow. In some 
cases, the wave Is more suitably modeled as a trochoid [9], l.e., 

?(_p) “ C(x,y) = Acos0— + sj + (3.13) 

and 

X = X - Asln(-7“ x+b). (3.14) 

o \A o / 

For a sinusoid the profile depends sinusoidally on x, while for a trochoid 

the profile depends sinusoidally on x^, which Is related to x by Eq. (3.14), 

a non-llnear transformation. The difference between a trochoid and a 

sinusoid can be seen from the following reasoning. For simplicity, let 

B = 0, so the wave crest occurs at x = x =0. We have 

o 


X > X 

o 

-A/2 < X < A/2 


X = X 

X = A/2, -A/2 

(3.15) 

o 

X < X 

-A < X < -Ml, A/2 < X < A . 
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Figure 3.5. Mean received pulse shape for reflection from an ocean wave 
with a height of 2 m (trough-crest). The laser altimeter is 
slightly off normal, with a nadir angle of 1“. 
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For points around the crest, > x, the decay of a trochoid from Its peak 
Is faster than a sinusoid, so Its peak Is sharper. For points near the 
trough, x^ < X and the rise from the trough Is slower, resulting In a 
smoother trough. In Fig. 3.6, we plot sinusoidal and trocholdal waves 
for a trough-crest wave height of 2 ni and a wavelength of 10 m. Indeed, 
the trochoid Is sharper around the crests and smoother around the troughs. 

Since trochoids have sharper crests and smoother troughs, we expect 
the reflections from the crests to be weaker and the reflections from the 
troughs to be stronger than In the case of sinusoids. This Is born out by 
the results shown In Fig. 3.7. For trocholdal waves, the reflections from 
the troughs are enhanced. These results should be compared with Fig. 3.3. 

Obviously, the ocean-surface profile Is not a perfect sinusoid or 
trochoid. The Gaussian model and the sinusoidal wave model can be viewed 
as the two extremes. The actual received pulse shapes are likely to lie 
somewhere between the broad smooth pulses, which result from the Gaussian 
model, and the highly structured shapes, which arise from the sinusoidal 
model. For the pressure-measurement application, a small laser beam 
divergence angle Is desired to minimize pulse broadening due to beam- 
curvature effects. Therefore, the received pulses are expected to exhibit 
fine scale structure due to strong reflections within the laser footprint. 
Later we will see that this fine structure is an important factor In 
achieving the picosecond (l.e., millimeter) timing accuracies which are 
required for measuring atmospheric pressure. 
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Figure 3.6. Comparison of sinusoidal and trochoidal waves. 
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Altitude = 400 km 

Transmitted Pulse Width = I cm (33ps) 
Wave Height = 2m (Trough- Crest) 

Beam Divergence = 100 firad 
Wavelength = 10 m 
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Figure 3.7. Mean received pulse shape from reflection from a trochoidal 

wave with a height of 2 m (trough-crest). The laser altimeter 
is at an altitude of 400 km. 
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4. TIME-RESOLVED SPECKLE 


4.1. Introduction 

For pressure measurements over the ocean, the ocean surface is the 
ranging target. Since the ocean surface is rough on the optical scale, the 
reflected laser pulses are corrupted by speckle. In this chapter, we study 
the statistical properties of speckle and its effects on the ranging 
accuracies of laser altimeters. The results not only are Important to the 
pressure measurement technique, but also have applications in general laser 
radar and ranging. 

Since the first discovery of the laser speckle phenomenon, there have 
been many studies treating its statistical properties. However, these stu- 
dies have been either for CW laser illumination [10], [11], or for a pulsed 
laser where the width of the received pulse is comparable to the correla- 
tion length of the speckle-induced fluctuations [12j. This latter case 
applies when the reflecting surface, or target, is rough on the optical 
scale but has a range spread which is much smaller than the laser pulse 
width. The speckle causes random fluctuations of the total received energy 
and is an Important noise phenomenon in target detection. 

With the advent of mode-locking and Q-switching technologies, laser 
pulses of a few picoseconds in durations are common today. The narrowness 
of the transmitted pulse promises higher accuracies in applications like 
remote sensing and ranging. In practice, most targets will have range 
spreads that far exceed the width of the laser pulse. Such would be the 
case for reflections of short laser pulses from airplanes, the ocean sur- 
face and the ground. The received pulse will be broadened to about twice 
the range spread of the target. The pulse shape is related to the geometry 
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of the target. In this case, the received pulse has a width longer than 
the correlation length of the speckle-induced fluctuations. As a con- 
sequence, speckle will cause random small-scale fluctuations within the 
received pulse which distort its shape. Similar phenomena can also occur 
for extended flat diffused targets, where the broadening of the received 
pulse is due to the wavefront curvature of the laser beam and for con- 
tinuously distributed targets in the atmosphere. A computer-generated 
example of a speckle distorted pulse is shown in Fig. 4.1. This phenomenon 
is referred to as time-resolved speckle. 

For applications like target identification and remote sensing of sea 
states [4], the waveform of the received pulse is used to characterize the 
target. Therefore, knowledge of the statistics of time-resolved speckle is 
important. In some other applications such as laser ranging and the 
pressure measurement technique considered in Chapter 2, estimation of the 
arrival time of the returned pulse in the presence of time-resolved speckle 
is the problem. 

In this chapter, we first derive the statistics of the detected 
signal. Then the problem of estimating the arrival time of laser pulse in 
the presence of time resolved speckle is considered. In Section 4.5, we 
study partially developed time-resolved speckle. A specific example is 
treated to illustrate the approach taken in the analysis. 

4.2. Statistics of the Detected Signal 

In this section, the statistical properties of the detected signal for 
laser pulses that have been reflected from a diffuse target are considered. 
The results are derived for a direct detection-type receiver system, which 
consists of a receiving telescope followed by a photodetector. 



TIME SCALE 


Figure 4.1. Example of time-resolved speckle. 
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Let the observation Interval be (0,T) and the observations be the pho- 
toelectron counts in N time bins each with width t . The observations will 

o 

be denoted by the photoelectron count vector k, and 


k *^2 * * * * * 


(4.1) 


where k^ is the photoelectron count within the i*"^ time bin. 

The statistics of k are related to the signal field a (r,t) received 

i s — 

CH 

over Che past seconds. Within Che i^ time bin, Che signal energy 
received by the system is 


Wi= / 


it 


‘It / 


(4.2) 


where we have assumed a simple aperture weighting function, which Cakes on 
the value 1 inside the aperture and 0 outside the aperture. is the 
receiver aperture area. For reflections from a diffuse target, a^(_r,t) 
follows circular complex Gaussian statistics [10]. 

We can also write W^ as an Infinite sum by performing a modal 
decomposition on the received field a^(£^, t). Using the Karhunen-Lodve 


expansion, we have 


a (r,t) = i(»^(r,t) , (i - 1 )t < t < it (4.3) 

s — mm— o — o 

m=l 

where form a complete orthonorraal set of basis functions, and are 

uncorrelated. Since a (r,t) has circular complex Gaussian statistics, Che 
coefficients are circular complex Gaussian variables, which implies 

that they are independent. 

By substituting Eq. (4.3) into Eq. (4.2) and making use of the ortho- 
normal properties of we have 
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“i = I 1^1^ • 

1=1 


(4.4) 


,th 


This shows that the signal energy received within the i time bin consists 
of contributions from infinitely many modes. The energy associated with 
each mode obeys the negative exponential distribution and is independent of 
the energies of the other modes. In the absence of photodetector satura- 
tion, the following relation holds: 


CO 



i 


n 


m 


(4.5) 


where n^ is the photoelectron count in the i^^ bin due to the m^^ mode. 

1 12 

The distribution of n when conditioned on lor Is Poisson, and since 

m ' m' 

12 1 
|o^j is negative exponentially distributed, n^ is a Bose'-Einsteln 

variable. The probability density function of kj, is then the infinite 

convolution of the probability density functions of n^, i.e., 


p(k^) = p ) h p(n2] * ... 


(4.6) 


Evaluation of Eq. (4.6) involves solving the eigen equation and carrying 
out -the infinite convolution. Simplifications are -possible with some 
approximations . 

For the l*"'^ time bin, the observation volume, which is the outer pro- 
duct of the receiver aperture area with the time interval, can be regarded 
as consisting of subvolumes or correlation cells [10], [13], with the 
energy density being approximately constant within each cell and statisti- 
cally independent of the energy densities of all other cells. The energy 
associated with each correlation cell is assumed to obey the negative expo- 
nential distribution, with the mean energy in each subvolume taken to 
be the same. This is equivalent to assuming the solutions of the eigen 
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equation in the Karhunen-Lo^ve expansion have only non-zero eigenvalues, 
each of the same value [14J, (15]. It turns out that for a simple aperture 
weighting function, this is a reasonable assumption [lOj. 

Under the above assumptions, the energy received during the 
i^^ time bin is gamma distributed, obeying the probability density function 




M, 


exp 

r w *1 
1 

-M ^ 

\ i 



L. ij 


(4.7) 


where is the mean value of the signal energy in the i^^ bln W^. 
Mathematically, this is nothing more than matching the actual probability 
density function of with a gamma density function. Since in practice, 
we will only be able to calculate the first and second moments of the 
intensity of the received signal, the most natural and simple way to choose 
the parameters of the approximating gamma density function is to match the 
mean and variance with the actual values. 

The expressions for the mean and variance of can be obtained by 
using Eq. (4.2) together with the properties of the circular complex 
Gaussian fields; the results are- [4], [6] 

_ iX y 

W = / ° dt / d"^r J (_r,t;£^,t) (4.8) 


and 


var(W^) 


It 

/ ° 

'(i-l)x^ 

o 


dt. 


/ "■ 
U-l)to 




1 ’^1 ’— 2 ’^ 2 ^ 


(4.9) 


where J is the mutual coherence function of the received field, 
a 

For reflections from diffuse targets, the mutual coherence function 
can be calculated using the Fresnel reflection formula together with the 
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properties of the circular complex Gaussian fields; the results are [4], [6] 
/ d^r J (r,t;r,t) = Q / d^p b,(p,z)jf(t - t|;) (4.10) 


and 


/ 


A 


R 


■*£1 / 

'Si 

2 

= -|- / d% b^(^,z) [f(t^ - <l>)l^if(t2 - 


(4.11) 


where Q is the expected received energy per pulse. The remaining variables 
were defined previously in Chapter 3. 

For a Gaussian shape transmitted pulse Intensity with rms width o^, we 
have the following simplified results. 


/ i '^a^^»^5I.,t) 


q F^(t) 


and 


(4.12) 


/ I 1-^a^-l • *^i 5-L2 ’ *^2 ^ I ^ 


qV^G(/ 2 (Jf,t^-t2) F^f— 


+ t. 


(4.13) 


where 


F^(t) 

F2(t) 



b2(^,z) G(o^,t - 


> 1 ') , 



(4.14) 

(4.15) 


G( a,t) 


a 


(4.16) 
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Fj(t) is the mean received pulse shape and F 2 (t) is the waveform of the 
variance of the received signal. Both are functions of the target. 

By substituting Eq. (4.12) into Eq. (4.8) and making the assumption 


that is much smaller than the width of F^(t), we have 

The variance of can be calculated from Eqs. (4.9) and (4.13) 

-2 it it It. + t.) 

var(W ) = -^ \ ° dt, / ° dt G(/2o ,t -t ) F [ = 

0 o 


(4.17) 


° dT^ G(/20j,T^) / 


2 

T, dT,F,(T,) 


+ / dTj G(/2aj,Tj) / 


-T 


T, dT, F(T,)J» 


(4.18) 


where we made the following changes of variables: 


T = t — t 
1 ^1 2 


and 


(4.19) 


= (t^ + t.p/1 . 


(4.20) 


Assuming the width of F 2 (t) is long compared with x^, and x^ is large com- 
pared with o^, we have 

var(Wj) -liFjdT,) z//''’ dt, G,( /Zd^.T, X - T,)| 




\ 


1 - e 





o o 


(4.21) 
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Finally, the parameter of the gamma density function can be obtained by 
matching the degree of freedom of the gamma density function with the 
signal-to-noise ratio. 


*^i “ Var(W.) “ F (it ) * 

1 Z o 

will be referred to as the number of 
^ th ^ 

1 time bln. 

If the laser footprint is uniform, 
to the laser footprint and the width of 


(4.22) 

speckle correlation cells in the 

or if the target is small compared 
the transmitted laser pulse is 


small compared to the range spread of the target, F^(t) and F 2 (t) given by 
Eqs. (4.14) and (4.15) are approximately equal and we have 


H K Fi(lTg)T^ . (4.23) 

Since ^ 1 is always smaller than K. Equation (4.23) is not a 

good approximation of Equation (4.22) in the leading and trailing edges of 
the received pulse where it could predict values of smaller than one. 
should always be greater than one. 

For the case when the receiver can not resolve the return waveform so 
only the energy of the received pulse is observed, the number of speckle 
correlation cells is simply K [121. For a target with large-range spread, 
only part of the target area is illuminated within one time bln and the 
area available for interference is smaller. As a consequence, the size of 
the speckle is larger and is smaller than K. 

Since is gamma dlstribtued with parameter the distribution of 
k^ can be shown to have a probability density function given by [12] 
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P(kj,) 


r(kj_ + 

r(k^ + i)r(M^) 



Here 


(4.24) 


k. =» nW /hv = <N> F,(it )t , (4.25) 

11 loo 

n Is Che quantum efficiency of Che phoCodeCecCor and <N> is Che expected 
number of phoCoelecCrons per pulse. When is an integer, Eq. (4.24) 
is the negative binomial density function. 

For most laser receivers Is longer than the transmitted pulse width. 
For this case, the photoelectron counts from different time bins are 
statistically independent. The joint density function of ^ is then 
N 

p(lc) = n p(k.) , (4.26) 

i=-l ^ 

This equation forms the basis for the analysis on receiver timing in 
Section 4.3. 

We will look at some specific target configurations. The first 
example is an extended flat diffuse target Chat covers Che entire laser 
footprint. For simplicity, we assume the target to have uniform reflec- 
tivity. This is a good model for reflections from Che ground, sands or 
walls of large buildings. In this case, the broadening of the received 
pulse is due to the curvature effects of the laser beam. 

Assuming a Gaussian laser cross section with a radius a^(exp(-l/2) 
point), we have 

■2a.l2 

For normal Incidence, the mean and variance of W^ can be calculated to give 
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(4.29) 


Equations (4.28) and (4.29) are plotted in Figs. 4.2 and 4.3 for three dif- 
ferent laser beam divergence angles. The number of speckle cells in the 
i^^ time bin is given by 

Kczt ..fczo, uU Tcz o, i r 1 

r o 


T - czo. it 

0^2 f o 

_ erfc = — 


(4.30) 


which is plotted in Fig. 4.4. The figure shows that starts increasing 
when the leading edge of the laser pulse illuminates the target; it remains 
constant after the trailing edge of the pulse arrives at the target. This 
can be seen from the formula of shown in Eq. (4.30). When the entire 
pulse illuminates the target, 1 is large and is approximately given by 




(4.31) 


which is constant with time. Equation (4.31) does not depend on so the 
curves of in Fig. 4.4 for three different beam divergence angles coin- 
cide with each other. 

For non-normal incidence, additional broadening due to the tilt of the 
target dominates the pulse shape. Neglecting the curvature effects, we 


find both the mean and variance of to be Gaussian in shape; they are 
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Figure 4.2. Mean received pulse shape for reflection from an infinite flat 

diffuse target. Normal incidence. Magnitude scale is relative. 
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Figure 4.4. Number of speckle correlation cells for reflection from an 
infinite flat diffuse target. Normal incidence. 
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plotted in Figs. 4.5 and 4.6. is plotted in Fig. 4.7. is constant 
in time. 

The next target we consider is a finite size flat plate. The 

geometry of the target is shown in Fig. 4.8. The face of the plate is not 

normal to the incident beam, instead, the plate is slanted with an angle 9 
and, therefore, has a non-zero range spread equal to the length of the 
plate times sin 9. This model can represent reflections from roofs, 
billboards or one facet of a vehicle. 

We shall assume the center of the plate is situated at the center of 
the laser footprint. For a rectangular plate of size £ m long by h m wide, 

the mean and variance of W^ can be calculated to give 




(4.32) 
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Figure 4.5. Mean received pulse shape for reflection from an infinite 

flat diffuse target. Non-normal Incidence. Magnitude scale 
is relative. 
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Figure 4.7. Number of speckle correlattoa cells for reflection from an 
Infinite flat diffuse target. Non-normal incidence. 
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with K 


K - 


When the plate Is small compared to the laser footprint, K Is simplified to 
A ah cos 0 

K = ■ j j . (4.35) 

X z 

The mean and variance of are plotted in Figs. 4.9 and 4.10 for three 
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Figure 4.9. Mean received pulse shape for reflection fron a finite size 
plate. Magnitude scale is relative. 
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Time Bin Width = 0.9 cm 
Target Distance = 300 m 
Plate Size = Im x 0.05m 
Receiver Aperture Area = 100 cm^ 
Laser Wavelength = 1064 nm 
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20. 60. 100. 140. 180. 
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Figure 4.10. Waveform of the variance for reflection from a finite size 
plate. Magnitude scale is relative. 
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different slant angles of the plate. The number of speckle cells is 
plotted in Fig. 4.11. The number of speckle cells starts increasing when 
the leading edge of the laser pulse illuminates the plate; it reaches a 
constant value when the trailing edge of the pulse arrives at the plate. 
The number of speckle cells starts decreasing when the leading edge of the 
pulse leaves the end of the plate. When the whole laser pulse Illuminates 
the plate, maintains a constant value given by 


Kct erf 
o 


«i = 


I cos 9 


2 a, 


2 

4/ir tan 0o^ erf^ 


i cos 9| 


2/2(3. 


/ir A ct a. erf — ^ 
R o 1 ' - 


2/2o, 


,22 „ h 

X z tan 0 erfly^ 


(4.36) 


For finite-size targets, the relative sizes of the target and laser 
footprint are important factors in determining the return pulse shape. For 
example, if the length of the plate is increased from 1 meter, which is the 
length of the plate used for Figs. 4.9, 4.10, and 4.11, to 3 meters, the 
corresponding plots become Figs. 4.12, 4.13, and 4.14, which are similar to 
the return waveforms from an infinite-size target. 

Another target we consider is a Gaussian random surface. This is a 
good model for reflection of laser pulses from the ocean surface. For 
negligible curvature effects, the expected received pulse shape is Gaussian 
[4]. If the transmitted laser pulse is short compared to the surface 
height variation, F^(t) and F^Ct) will be approximately equal, and 


NUMBER OF SPECKLE CELLS 


48 


TIME BIN NUMBER 


Figure 4.11. Number of speckle correlation cells for reflection from a 
finite size plate. 
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Time Bin Width = 2.7cm 
Target Distance = 300 m 
Plate Size = 3m x 0.05m 
Receiver Aperture Area = 100 cm^ 
Laser Wavelength = 1046 nm 



0 . 40 . 80 . 120 . 160 . 200 . 

20 . 60 . 100 . 140 . 180 . 

TIME BIN NUMBER 

Figure 4.12. Mean received pulse shape for reflection from a finite size 
plate. Magnitude scale is relative. 
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Target Distance = 300m 
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Figure 4.13. Waveform of the variance for reflection from a finite size 
target. Magnitude scale is relative. 
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(4,37) 


where is the rras height of the surface profile. In this case, is 
roughly proportional to the pulse shape. 

In summary, the received signal in the 1^^ time bln is gamma distri- 
buted with parameter M^, and when is an Integer, the photoelectron count 
is negative binomial distributed. is the slgnal-to-nolse ratio of the 
received energy and is referred to as the number of speckle correlation 
cells. From the examples considered, we find that the mean received wave- 
form is closely related to the geometry of the target. The statistics of 
time-resolved speckle also depend on target geometries. 

4.3, Estimation of Arrival Time 

In laser ranging and altimetry applications, the arrival time of the 
returned pulse, denoted by t^, is of interest. In this case, Eq. (4.26) 
can be written explicitly as 

kl( T^)\-«i( Tj) 

wj 

(4,38) 




N 

n 

i=l 


r(k^ + m^(t^)) 
r(k^ + i) r(M^(r^)) 


1 + 


MjT^\-ki 

ki(x^)/ 


The Maximum Likelihood (ML) estimate of the arrival time is the value of 

T, that maximizes the likelihood function p(tt|T ), i.e., 
d d 


T, = arg max lp(klx.)] 

ML T. 

d 

= arg max [ £n p(k|x^)l , (4.39) 

^d 

where arg max[F(x)] denotes the argument x that maximizes F(x). By substl- 

X 

tuting Eqs. (4.38), (4.25) and (4.22) into Eq. (4.39), we have 
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T, = arg max 

T 

a 

N 


I tor(k, + M (t,)) - I £nr(k, + 1) - I tor(M,(T,)) 
i=l i 1 d i iil id 


KF,(iT - T.) 
I o a 


N 


<N>F (it - T )' 
z o a 


- k^ Zn^l + -f 

(4.40) 

Assuming that the received pulse always stays Inside the observation inter- 
val for all values of so that there are no end effects, we obtain 


arg max 


ML 


oax j; (kj+Mj^(Tj) - y) to(k^+M^(T^)) - I - jjin(M^(T^)) 

d ^ ^ 


N 


/ KF (It - r ) \ N / <N>F (iT - T. 

- + <H>F^(lt^ - T^)j - - T^> 


(4. 


41) 


where we have made use of the following relation 


inr(x) 


- (* - t) 


tox - X + Y to2ir 


X > 1 


(4.42) 


Assuming the summations of M (t ) over i to be independent of t , we have 

id d 


arg max 


ML 


1 

IX I (k^ + MjCt.) - y) ^(k. + M (t.)) 

TdU=l 


- I k to 1 + 

i=l ^ ' 


- -^d> 

<N>F2(i^o - -^d>> 


(4.43) 


When the footprint is approximately uniform over the target area and the 
width of the laser pulse is short compared to the range spread of the 
target, Fj^ and F^ are approximately equal. Equation (4.43) reduces to 


T, = arg max 
‘^ML T 


d ^ 


Jj W - 7 ) + V’d» 


= arg max(H(T.)] 
^d ' 


(4.44) 
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The performance of t, will be analyzed based on a procedure similar 

“ml 

to that used by Bar-David (161. 

If H(Tj) is differentiable, we have 
d 


d 




H( T. ) = 0 . 


(4.45) 


If the actual arrival time of the pulse is denoted as then H( t^) can 

be expanded using the Taylor series expansion centered around as 


, " ”^do^ 

a( T^) - H( + H( ^d " ■'do^ ^do^ 2l • * 


(4.46) 


Differentiating Eq. (4.46) with respect to and using Eq. (4.45) yield 


H(t ) - 0 - H(t ) + (tj - H(tj ) + ... , 

ML ML 


(4.47) 


from which it follows that the error of the estimator is 


e = T , - T, 

‘‘ml 


« "do> 


(4.48) 


provided H( t, ) * 0. This constraint excludes strictly rectangular pulses 
do 

from this analysis. 

The bias and variance of the estimator can be written as 


E[£] = -E[H(T^^)/a(T^^)l 


(4.49) 


and 


var[ e] = E< kr 


fr«‘ 'do>i 


1 

o 

yy 

i 


- E^(e] , 


(4.50) 


where the expectations are over the observation vector k. 


Assuming )] » var[H( t, )], such that H( x , ) can 

do do do 

by its expected value, we have 


E[e] 


E{H*(x^^)} 


and 


var[ e] 




ElH(Xdo)] is given by 
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N - 


B[H<Xj^)I . I E«n(k^ -X 1 -k-nr 

i=l I i=l i i 


N 




To proceed further, we note that 


Ak. 


£n(k, +M.(t, )) = to(k. +M,(x, )) + £n 1 += ^ 

i 1 do i i do \ k + M.( x, ), 

1 1 do 




= ta(k. + M. (x.„)) + = 


. T i-i. V ^ = ; T > 

X X uO If + M r T ^ 

*^i ”i’^ Mo^ 


where we have separated k^ into a deterministic part k^ and a 
random part Ak^. The validity of Eq. (4,54) requires 


be replaced 


(4.51) 


(4.52) 



(4.53) 


(4.54) 


zero mean 
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which is satisfied when the received signal level is high or when the 
number of speckle cells is large. 

Substituting Eq. (4.54) into Eq. (4.53) yields 
N 

E[H(t, )] 5 K T I F (it - T ) in F (it - T ) 
do o 1 o do . 1 o do 


s Kt^ / F^ in Fj = 0 ; 


(4.56) 


therefore, the estltaator is unbiased. The variance of H( t, ) is 

aO 


var[H(Tdo)l =■ I 


N 

I 

' 1=1 


«i^^do> 


‘^1 '‘i^^do^ 


var( Ak^) 


(4.57) 


From the probability density function of k^ shown in Eq. (4.24), we have 
•- "^1 

var[Ak^] “ ’‘i TT* ‘ (^.58) 

By substituting Eqs. (4.58), (4.25) and (4.23) into Eq. (4.57), we obtain 


N F? 




(4.59) 


E[H(Tdo)l is given by 


N 


E[H( 






N .. 


N ^do> 


y M. (t, ) in(k. + M.(t, )) + y = 

1-1 ‘ i ‘ ■io iiikj(ta„) + "i<W 


“i N F? 


_ I IL, 

Un> K/ F^ o ’ 


(4.60) 
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where use has been made of Eq. (A. 54). Finally, by substituting these 
results into Eq. (4.50), we have for the ML estimator 


var[ 



1 


N 

I 


i=l 



T 

O 


(4.61) 


As a comparison, the ML estimator derived for the case of negligible 
speckle noise is [13], [16], 


arg max 


£n 


‘d 


N 


I £n k^( T^) 


i=l 


(4.62) 


It will be interesting to see how this estimator performs in the presence 
of speckle. Employing the same procedure used above, we find it to be 
unbiased and the variance of its error is given by 


var[ e] 
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2 

^ 1 
1 i=l 


1 2 


<N> 


T 

I T~ 

i=l 1 


■ (J, S’ 


(4.63) 


which reduces to Eq. (4.61) when and F^ are equal. This Indicates that, 
under the assumptions made in analyzing the two estimators, l.e., shot 
noise and speckle noise are not severe, their performances are 
indistinguishable. 

It can be seen from Eq. (4.63) that the timing error is smaller for 
waveforms with larger first derivative. The effects of signal bandwidth on 
the timing performance will be further discussed in Chapter 5. 

We shall use Eq. (4.63) to calculate the timing performances for the 
target geometries considered in Section 4.2. For normal incidence from an 
infinite flat target, the mean received waveform corresponds to Fig. 4.2. 
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The ranging accuracy is calculated and plotted in Fig. 4.15 versus the 
average number of photoelectrons. From the figure, we see that in the 
limit of high signal energy all three curves have the same accuracy limit 
set by speckle. This is expected because they have the same number of 
speckle cells as shown in Fig. 4.4. For non-normal incidence, the mean 
waveform was shown previously in Fig. 4.5 and the ranging accuracy is 
plotted in Fig. 4.16. For both normal and non-normal incidence cases, 
geometries that give rise to sharper received pulses result in better 
timing accuracies. 

For the finite-size plate target considered in Section 4.2, the return 
waveform was found to depend on the relative sizes of the plate and the 
laser footprint. Therefore, we expect the ranging performance will also 
depend on the size of the plate. For a plate length of 1 meter, the mean 
waveform corresponds to Fig. 4.9; the ranging performance is calculated and 
plotted in Figure 4.17. It is interesting to note that, although for a 
slant angle of 60“ , the return pulse is broader than that of a slant angle 
of 35®, the ranging performance is actually better for the same signal 
level. This can be explained by examining the slopes of the return 
waveforms shown in Fig. 4.18. For the 60® case, the ends of the plate are 
located closer to the center of the footprint, so when the laser pulse 
first hits the front end of the plate, the Increase of the reflected energy 
is larger. Therefore, the slope of the leading edge of the reflected pulse 
is larger. For the 35" case, the ends of the plate are farther away from 
the center of the footprint, so the slope of the leading edge of the 
reflected pulse is smaller. From Eq. (4.63), we see that waveforms with 
larger slopes or bandwldths result in better timing accuracy. We need to 
keep in mind that in Fig. 4.16 the performance is plotted versus the 
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Figure 4.15. Timing performance for reflection from an infinite flat 
diffuse target calculated using Eq. (4.63). Normal 
incidence. Mean received pulse shape corresponds to that 
in Figure 4.2. 
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Figure 4.17. Timing performance for reflection from a finite size 

plate calculated using Eq. (4.63). Mean received pulse 
shape corresponds to that in Figure 4.9. 
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average number of photoelectrons. For a fixed laser output energy and 
fixed plate size, the average number of photoelectrons will decrease as the 
slant angle increases. In practice, this has to be taken into account in 
comparing the timing performance. 

For a 3 m long plate, the return waveform is shown in Fig. 4.12, and 
the slope of the waveform is plotted in Fig. 4.19. The ranging perfor- 
mance is shown in Fig. 4.20. In this case, a 45“ slant angle results in 
the worst timing performance. Again Figs. 4.21 and 4.22 clearly show that 
waveforms with larger derivatives permit better timing accuracy. 

For a Gaussian received pulse shape with rms width Eq. (4.63) can 
be evaluated to give 

varlc! “ 7)4 ’ (4.64) 


where we have approximated the sununation by an integral. If <N> » K, then 
the timing accuracy is limited by speckle; on the other hand, if K » <N>, 
the timing accuracy of the receiver is shot noise limited. 

A suboptlmal estimator that could be simpler to Implement is 


arg max 


d*- 


ill “1 ''i' '1' 


(4.65) 


which is the discrete time analog of the correlation receiver or matched 
filter. Again, the estimator can be analyzed by Bar-David's [16] method. 
It is easy to show that the estimator is unbiased, and its error has a 
variance given by 



(4.66) 
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.19. Slope of the mean pulse shape shown in Figure 4.12. Magnitude 
scale is relative. 
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When both and ?2 are equal to a Gaussian pulse shape with rms width o^, 
we have 

(<k * i) 

Comparing Eq. (4.64) with Eq. (4.67), we find that has a rms 

c 

timing error which is about 25% greater. In critical timing applications, 

it may be necessary to use the better estimator t. . 

“in 

From the above results, we can see the effects of speckle on the 
receiver timing performance. When speckle is present, the mean-square 
timing error is increased by a factor of <N>/K. If K is much larger than 
<N>, the effects of speckle on timing are negligible. But if <N> is 
comparable to or even greater than K, speckle will degrade the timing 
performance significantly. 

4.4. Computer Simulation 

Since Eqs. (4.63) and (4.66) were derived under the high signal-to- 
noise assumption, they are not accurate for low signal levels; other 
methods must be used to determine the estimators' performances for weak 
signal returns. A computer simulation of maximum likelihood estimators has 
been developed by J. B. Abshire and J. F. McGarry at NASA Goddard Space 
Flight Center (GSFC) for this purpose. This simulator was an extension of 
earlier work [17] used to evaluate the low-level timing performance of 
estimators with shot noise only. A summary of simulation results was 
recently presented [18]. 

The simulation is performed by dividing the observation interval into 
512 time bins and using the computer to generate the photoelectron counts 
in each time bln. The probability density function of the photoelectron 
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count In the time bin pCk^^) was given previously in Eq. (4.24). 
integer value of M^, p(kj^) becomes 



For an 


(4.68) 


which is the negative binomial or Pascal probability density function. 

The negative binomial distribution originates from a sequence of inde- 
pendent Bernoulli trials. In such trials, each trial has a fixed proba- 
bility of success, p, and fixed probability of failure, q. The probability 
of k failures prior to the r*"^ success 1 19] is given by 


p(k) 



r 

P 



(4.69) 


This probability density is called negative binomial. The derivation of 
this formula is straightforward after observing that the r^^ success 
occurs on the (r + k)*”^ trial if and only if (r - 1) successes and k 
failures occurred in the first (r + k - 1) trials and they are followed by 
a single success. 

Comparing Eq. (4.68) and (4.69), Abshire and McGarry found that the 
probability of observing k^ photoelectrons in the 1^^ time bin was equal 

t* h 

to the probability of k^ failures prior to the success in a sequence 

of Bernoulli trials with probability of success equal to + k^). 

Using this analogy, they generated k^^ by conducting a sequence of Bernoulli 

trials and counting the number of failures which occurred before the 
cH 

success occurs. The advantage of this approach is that the Bernoulli 
trial can be Implemented easily and efficiently on a computer. The computer 
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generates uniform random numbers between 0 and 1. Whenever Che generated 
number Is less than a threshold, the trial is classified as a success. The 
threshold is chosen to be Che same as the probability of success, which is 
Mf/ (Mf + k^) for the 1^*^ time bin. 

In general situations, may not be an integer and Eq. (4.24) does 
not reduce Co the negative binomial distribution. In these cases the simu- 
lator rounds to the nearest integer value. This quantization of only 
introduces small errors if is large. 

The simulator was run for a raised cosine waveform, which can be 
expressed as 

F(t) - 1 (l + cos ^ t) ; -jltly . (4.70) 

is calculated using Eq. (4.23) and rounded to the nearest integer. If 
falls below 1, it will be replaced by 1. The results of the simulator 
are compared with Che theory in Fig. 4.21. The agreement between Che 
theory and simulation is better for high-signal level and for large K. 

This is expected, because in deriving the theory, we have assumed a high 
signal-to-noise ratio, so that the approximation that leads to Eq. (4.63) 
is valid. 

For K equal to 5040, the transition occurs at signal energy about 1000 
photoelectrons, where the simulation and the theory come into agreement. 

4c this signal level, the average number of phoCoelectrons in each time bin 
is only about 2. In practice, the signal level is usually higher and Che 
theory is expected to predict the timing performance accurately. 

The simulator has two limitations. The first is that in its present 
form it can only work with Integer values of 


In principle, we can 
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generate the counts for a non-integer by first generating a gamma 
distributed variable, then use it as the mean function in a Poisson random 
variable generator, and generate the desired variable. Although this 
approach requires more computer time, it gives more accurate results for 
low values of This extension of the simulator is presently underway at 

GSFC. The second limitation is that the timing resolution is limited by 
the finite number of time bins (presently 512) in the simulator. 

Therefore, the arrival time of each simulated waveform can only be an 
Integer multiple of the width of a time bin. 

The error introduced by the timing quantization can be calculated by 
using the fact that it is uniformly distributed over one time bin. 

Therefore it has a standard deviation of 0.29 ( 1/ /12) time bin. This is 
equivalent to an error of 5.6 x 10 ^(0.29/512) when normalized by the pulse 
width. The quantization noise is not Important for the results shown in 
Fig. 4.21, but it is the limiting factor when the timing error is small. 

Therefore, if the simulator is to be used for waveforms that result in a 

normalized timing accuracy better than 0.0005, more time bins are required. 

Unfortunately, an increase in the number of time bins is usually accom- 

panied by a quadratic Increase in the computer execution time. This execu- 
tion time is the limiting factor with the present simulator. 

Since the computer simulation results and the theory are In good 
agreement for high signal levels, this verifies the theory’s validity in 
that region. The simulator establishes the ML algorithm performance at low 
signal levels. 

4.5. Partially Developed Speckle 

In previous treatments of the time resolved speckle, we have assumed a 
fully developed speckle, i.e., the speckle follows the circular complex 
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Gaussian statistics, which Is generally true for diffuse targets or for 
targets with many reflection centers. For some targets, the speckle Is 
only partially developed, and the distribution of the speckle Is not 
Gaussian. An example of the partially developed speckle Is laser reflec- 
tions from the LAGEOS. 

The Laser Geodynaralc Satellite (LAGEOS) was launched in May 1976 and 
was designed as a passive long-lived target with a well-defined orbit [20]. 
As such. It functions as a reference point In the Inertial space and by 
ranging to It, sets of ground-based laser systems may recover their Indivi- 
dual geometry, their positions with respect to the earth's center of mass, 
or their positions with respect to an Inertial reference. In order to 
enhance Its reflectivity as a laser target, the satellite Is covered with 
optical cube corners which retroreflect any Incident optical signal. There 
are a total of 426 cube corner reflectors. The LAGEOS has a diameter of 
60 cm and Is orbiting at an altitude of 5000 km. Simple calculations show 
that, even for a reasonably large receiving telescope, say 1 meter In 
diameter, the Individual retroref lector on the LAGEOS can not be resolved. 

During the prelaunch testing at the NASA-GSFC In December 1975 and 
January 1976, there had been concerns about the pulse-shape fluctuations 
due to the coherency effects [20] . This Is due to the fact that there Is a 
random phase associated with the optical field reflected from each 
retroref lector . When the reflected fields overlap In time at the receiver. 
Interference causes the returned power to vary randomly. However, 
receivers with enough resolution to experimentally observe this pulse-shape 
fluctuation or time-resolved speckle were not available at the time of the 
prelaunch testing. The purpose of the work here Is to provide a theoretical 
analysis of the statistics of this partially developed time-resolved speckle. 


72 


Using the Fresnel diffraction formula, the optical field on the plane 
of the receiving telescope is 


aj^(jr,z,t) = c^ exp 


2 ~ 
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-Ik r^ 
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1=1 

o^i — 

2z 

L z J 


exp[i(|)^I 


(4.71) 


where 


c = (T /X z) exp(-i2k z] 
o a o o 


(4.72) 


and 


A^(t) 
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z) f 




exp 




(4.73) 


where a^(£^, z) is the complex amplitude cross section of the laser foot- 
print, is the Intensity transmittance of the atmosphere, is the 

wave number, and X is the wavelength of the laser. N is the total number 
o r 

of retroref lectors , 8^^ is the effective reflecting area of the 

retroref lector, which is related to the angle of the incident laser 

beam and the orientation of the retroreflector . In practice, most of the 

8^'s are zero due to the shadowing effects. is the random phase angle 

associated with the reflection from the i^^ reflector; we assume to be 

statistically Independent of each other and uniformly distributed between 0 

til 

and 2tt. ^ is the transverse coordinate of the i reflector and is its 
displacement relative to the reference plane. 

For a direct detection receiver, the received power is given by 


N N 



(4.74) 
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where *^(_p) is the two-dimensional Fourier transform of W(r), the aperture 
function. 

Since the individual reflectors on the LAGEOS cannot be spatially 
resolved, Eq. (4.74) can be simplified to 


P(t) . Icj2 A^ 


N 


I A. (t) exp[i<j> ] 
i=l 


(4.75) 


where we have made use of the following relationship, 

/k 




(4.76) 


Since are random, P(t) is also random. This fluctuation of the 

received power is due to the coherent additions of fields with random pha- 
ses. Because the number of components is not large enough to warrant the 
use of the central limit theorem, the statistics are different from that of 
a fully developed speckle. The probability density function of P(t) is 
derived in the following. 

Let the dummy variables x^, y^^, X and Y be defined by Che following 
equations , 


Xj^(t-) = A^(t) cos , 

yj^(t) = (ICgl^ A^(t) sin , 

N 

I X 
i=l ^ 


(4.77) 

(4.78) 

(4.79) 


and 


^r 

I Yj 
1=1 ^ 


(4.80) 
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Without loss of generality, we shall assume all Aj^'s(t) to be real. 

are uniformly distributed, the probability density function of 
be obtained by a simple change of variable. The results are 


^*i^ ”‘7 \ ; l*il <- l<^ol 


otherwise. 


the corresponding characteristic function is 


4» (m) “ / dx p (x ) exp[-iaK ] = J ( [c [Ay^A (t)ui) 

- ^ i ~i i o'o'Ri. 


til 

where J (a) is the Bessel function of the first kind, zero order, 
o 

Using the fact that <^^'s are independent, we have 
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Similarly, we have 
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The joint characteristic function of X and Y Is defined as 


^ (u,v) = / dX / dY p (X,Y) exp[-i(uX + vY) ] . 

X,Y — -« X,Y 

Making the following changes of variables 


~2 ~2 1/2 
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9 = arctan(Y/X) , 
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C = (u + V ) 
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x^ can 
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(4.81) 
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and 
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3 = arctan(v/u) , (4.89) 


we have 

«» 2ir 

<i> (u,v) = / dr / d0 p 0 ( 1 :, 0) exp[-irc cos(3 -0)] . (4.90) 

XY o o 

Since X and Y have Che same distribution, It Is reasonable Co assume chat r 
and 0 are circularly symmetric, so that r Is independent of 0 and 0 is uni- 
formly distributed between 0 and 2ir. Under this assumption, we have 


(j) (u,v) = / dr -p (r) J (?r) = / . dr p (r) J (r/i-Tv ) . (4.91) 

vv" ° \ r o 

XY o o 

Letting v equal 0 in Eq. (4.91), we obtain 


<|i (u,0) » (|i (u) = / dr p (r) J (ur) 

XY X" o ° 


= / drr[p (r)/rjj (ur) . (4.92) 

o 

(u) is recognized to be the Hankel transform of p (r)/r. Using the 
X ^ 

inverse Hankel transform, we get 

00 

.p (r) = r / duu 4 (u) J (ur) ., r 2. 0 
o X 


= 0; r <0 . (4.93) 

Finally, since P(t) is the square of r(t), we have 
00 

pp(P) “ / duu 4) (u) J (uP^ ^) , P >_ 0 

o X 

=0 ; P < 0 , (4.94) 

which is equivalent to Eq. (63) of BarakaC [21]. 
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Using Eq. (4.83), the probability density function of optical power 
can be written explicitly as 

r N 1 


pp(P) = (1/2) / du 
o 




u J^(uP^^^) ; P 2. 0 


0 ; P < 0. 


(4.95) 


In principle, Pp can be calculated for all t and for all possible ranging 
geometries, although the computation can be complicated. In practice, we 
are Interested In the statistical properties of the detected signal, and 
usually knowing the first- and second-order moments of P(t) is sufficient. 
In this case. It Is not necessary to evaluate Eq. (4.95). 

The mean power can be obtained from Eq. (4.75) directly: 

N 

E(P(t)] = Ic A_ I lA(t)l^ . (4.96) 

o I 


The autocovariance function of P(t) Is given by 



where we have made use of the fact that the retroref lectors cannot be 
spatially resolved by the telescope. 

Let the Impulse response of the photodetector be h(t) and its quantum 
efficiency by n; the mean and covariance functions of the output S(t) are 


given by 
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where 

p[ 2q 

t|) -2z/c-~ + — 2^ . (4.102) 

cz c 

p Is the coordinate of the center of the satellite. The first term In 
Eq. (4.101) Is due to shot noise while the second term Is due to speckle. 

We see that shot noise has a correlation length of about Z/2 and 
speckle has a correlation length of about . Also, we find 

that for the second term In Eq. (4.101) the terms Inside the double sum- 
mations are significant only when and satisfy 

< 2/2 Of (4,103) 


or 



5 ^) < 2/2 


(4.104) 


2 

Since p^/cz, the optical beam curvature term. Is very small, we have 
equivalently 

- 5j) < Op . (4.105) 

This shows that only the reflections from reflectors that are separated 


vertically less than the width of the transmitted pulse will interfere with 
each other. For a transmitted pulse width of 24 psec, only reflectors that 
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are separated less than 2 cm will cause interference, but for a pulse width 
of 2 nsec, reflections from reflectors as far as 42 cm apart will Interfere. 
Therefore, to minimize speckle, a shorter laser pulse should be used. In 
this case, only when two or more retroref lectors happen to lie on the same 
plane will there be intereference. 

We plot the waveforms for two possible LAGEOS orientations and the 
results are shown in Figs. 4.22 and 4.23. In the figure, the rras widths of 
the transmitted pulse and receiver impulse response are both equal to 1 cm. 
The time and the magnitude scales in the figures are relative. From the 
figures we see that, while the variance due to shot noise is proportional 
to the signal, the variance due to speckle may not be. For example, in 
Fig. 4.23, the leading edge of the return pulse is caused by the reflection 
from only the first reflector; no interference occurs so there is no 
speckle. In Fig. 4.24, we show the results for the same geometry as in 
Fig. 4.23, but both the widths of the transmitted pulse and receiver 
Impulse response are doubled. In this case, the long transmitted laser 
pulse smooths out the received pulse, and the return from the individual 
reflector cannot be temporally resolved. 

4.6. Summary 

In this chapter, we studied the statistical properties of time-resolved 
speckle and the effects of time-resolved speckle on the timing performance 
of the receiver. The results Indicate that, in general, the presence of 
time-resolved speckle causes the mean square timing error to increase 
roughly by a factor of <N>/K. If K is smaller than <N>, the timing 
accuracy is limited by speckle. 

Since the theory is derived under the high slgnal-to-noise ratio 
assumption, the timing performance of the receiver in the low signal level 
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Laser 


Figure 4.22a. Schematic diagram of the laser and the retroref lectors on 

LAGEOS. Relative position in units of transmitted pulse width. 
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Laser 


Figure 4.23a. Schematic diagram of the laser and the retroref lectors on 

LAGEOS. Relative position in units of transmitted pulse width. 
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Laser 


Figure 4.24a. Schematic diagram of the laser and the retroref lector on 

LAGEOS. Relative position In units of transmitted pulse width. 
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case can be more accurately obtained by a computer simulation. The results 
from the theory and simulation come into good agreement when the slgnal-to- 
nolse ratio is high. 

There are targets that give rise to partially developed speckle. The 
reflections of laser light from a specific target, the LAGEOS, are investi- 
gated to illustrate the approach taken in analyzing partially developed 
time-resolved speckle. 
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5. ESTIMATION OF DIFFERENTIAL ARRIVAL TIME 


5.1. Introduction 

To infer pressure, we need to measure the difference in the round-trip 
propagation times between two pulses which were transmitted simultaneously 
at two different wavelengths. In this chapter, we first evaluate and com- 
pare the performances of various timing algorithms. The relation between 
the timing accuracy and signal bandwidth is investigated in Section 5.3. 

In Section S.4, timing algorithms are compared for simulated ocean 
reflected pulse shapes. In Section 5.5, the transmitter effects are con- 
sidered. The effects of the sampling process are studied in Section 5.6. 

5.2. Timing Algorithms 

The Maximum Likelihood (ML) estimator for the differential delay can 
be expressed as 


■ 'iV 


^2 “ 

ML^ ML2 

arg max(p(Sj(t),S 2 (t) i T^)) - arg max[p(S^(t) .S^Ct) j x^) J , 


(5.1) 


where Xj^ and x^ are the arrival times of the laser pulses at the two wave- 
lengths. In Eq. (5.1) S^(t) and S 2 (t) are the received signals at wave- 
lengths and A^, respectively. '^ 1^2 ^ 2 f ^2 ^ estimates of 

x^ and X 2 when both Sj^ and S 2 are observed. 

Since the two wavelengths are widely separated, the speckle-induced 
fluctuations at the two channels are uncorrelated. Also, since the shot 
noises of the two channels are Independent, Eq. (5.1) can be simplified to 
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arg max[p(S 2 (t) | T 2 ) 1 
”^2 


- arg max[p(S^(t) | t^) ] 


= T, - T, , (5.2) 

ML ML 

where ^ 2 ^^ and are the ML estimates of individual channels. Equation 

(5.2) shows that the ML estimate of the differential delay is simply the 
difference of the single channel ML estimates of the arrival times of the 
two wavelengths. 

The single channel ML estimate of the arrival time when speckle is not 
severe was given previously in Chapter 4; the integral version of the esti- 
mator [13] is written here for convenience 

A 

T. = arg max 
nL 

where is the expected value of Unfortunately, Implementation of the 

ML estimator requires prior knowledge of the expected received pulse shape. 
Because of the random and dynamic nature of the ocean surface, the received 
pulse shapes will change randomly as the surface profile within the 
footprint changes. Since it is not possible to predict the expected pulse 
shapes a priori, a suboptlmal estimator that does not require knowledge of 
the pulse shape is needed. 

If the laser beams at the two wavelengths are aligned so that their 
footprints overlap, the reflected pulse shapes will be almost identical. 

We can write and S 2 as 

T^(t) = <N^> F(t) (5.4) 


" «0 

/ dt S^(t) in Sj^(t + T^) , 


(5.3) 


and 
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S2<t) = <N^> F(t - 


(5.5) 


where F(t) Is the normalized mean received waveform. One of Che simplest 
techniques for estimating Che differential propagation time is to calculate 
the peak of Che correlation function of Che two received pulses 


T = arg max 

T L-« 


/ dt S^(t) S^Ct + t) 


arg max » 


(5.6) 


where Rj ^2 the cross correlation function. 

Assuming R^ 2 ( t) to be differentiable, we have 


Rl2(ti2) 


(5.7) 


By expanding Rj 2 (t) in terms of the Taylor series expansion centered around 


Che actual differential delay Tj 2 » we can write [16] 


^12 ■’ ”^12 


Ri2( Tj2) 


(5.8) 


where e is the error of the estimator. 

If Che shot noise and speckle are not severe so that 
e{Ri 2 (Ti 2 )} » Var (r^ 2^ ■^12^ Square Error (MSE) 

are given by 


Bias = E{e} 


E{R^,(t^2>> 

E{Ri2(Ti2)} 


(5.9.) 


and 


MSE = E{e } = 


E{Rj2(Tj2>> 
E{Ri2(t^2)}^ 


(5.10) 
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The numerator on the right-hand side of Eq. (5.9) can be written 
explicitly as 

CO 

3 


E{Ri2(Ti2)}= / dt S2<t + t)| 


T=T, 


12 


<Nj><U 2> / dt F(t) F(t) . 


(5.11) 


Equation (5.11) is zero if the entire return pulse always stays inside the 
observation interval. Therefore, the correlation algorithm is unbiased. 
The MSE given by Eq. (5.10) can be evaluated as follows: 


E{R"2<Ti 2)} = Ei 




— 00 —00 


"r'^12 




^2 " '^121 


dt^ / dt2 Rg^(tpt2) ^82^^! '^2^| 


— CO —00 


T a T =» T 
1 ^2 ^12 


and 


m 2 

E{R(Tj 2^^ '“ / dt S'^(t) -^-■S2.(t + t)I 

-a" 3^ X 


T=T 


12 


where 


(5.12) 


(5.13) 


Rg (tpt2) = E{S^(tj) S^(t2)} = Cg (tj,t2) + S^(t^) S^(t2) 
By substituting Eqs. (5.12) and (5.13) into Eq. (5.10), we have 


(5.14) 


-<o -00 j 12 2 


MSE = 


-r = T = T 
1 2 12 


/ dt S^(t) S2(t + t^2> 


2 




The first terra on the right-hand side of Eq. (5.15) equals the square 
of the bias, which was shown previously to be zero. The second and third 
terras are the first-order terras. They represent the errors from corre- 
lating the received signals with their corresponding mean waveforms. The 
last term is the second-order term. It is the additional error from corre- 
lating two fluctuating pulses. If the mean pulse shape were known in 
advance and each received pulse were correlated with its mean waveform, 
then the sum of the second and third terms would be the total error. Since 
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we do not know the mean pulse shape, we correlate the two received pulses 
with each other. In this case both pulses are distorted by speckle and 
shot noise and correlating two distorted pulses gives rise to the addi- 
tional second-order error term. 

For a target that gives rise to a Gaussian shape mean received wave- 
form, Eq. (5.15) can be evaluated in closed form. While the detailed 
calculations are shown in Appendix A, the results are 


MSE 




<N^> 


2/2(a^ + 


+ 


<^f ^ 


2.3 

Oj.) 


(2<^^ 




2.3/2 , 2 3 2^3 2,3/2 

Op) ( «h + T + T Op) 



2/2 


2 / 2 

0, 

T h 


+ + 


4 )' 

T 


<24 


+ 2 a- 


^)3/2 


<4 


■^T 


a- + 






2.3 


<Ni>^ 


<N2>Kj 


2/T (a^ + 2ap(.^ + a^ + a^) 

(2a^ + <-%+Y 4'^ 


<N^XN2> 


0/ 2 . 2,. 2 ^ 2 . 2,3 

2(Op + OpXoj + <^ + oj) 

3/ 2.-2.- 2,3/2 


,2,2 

207 ( 0 ^ 


4 


2,3 

a^) 


K 1 K 2 


^4 


2 . 3/2 

0f) 


^4 


4 


20 ^'/' 

T 


(5.16) 


where a^ and are the rms widths of the transmitted pulse intensity and 
receiver impulse response, respectively, and 0 ^ is Che rms time spread of 


Che target. 
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For the case of reflection from a very small or point target, 0 ;^, = 0, 
Eq. (5.16) reduces to 


MSE 


(<N^> <N2>) 


2/2 aj((^ + 


1^1. 2.ZH .2.3 23/2 

(2aj^+Of) 


(<N^>K2 <N2>Kj) 


2 ^ oj(o^ + oj)^ 


2 ^ 2,3/2 , 2 3 2,3/2 

(20h+Of) (Oh-^y^f^ 


-2.2^ 2,3 

<N,XN,> 3, 2 ^ . 2,3/2 

1 2 Oj^(Oh + 2<^f) 


(5.17) 


Here we see that speckle does not contribute to the first-order timing-error 
terms. In this case, the target does not cause any broadening in the 
return pulses. Therefore, the speckle-induced fluctuations have a correla- 
tion length as wide as the received pulse. The effect of speckle is only to 
cause the amplitude of the received pulse to fluctuate and not distort the 
received pulse shape. Since the speckle only affects the amplitude of the 
received pulse, it will not degrade the correlation. However, shot noise 
still causes distortion of the received waveforms because it has a correla- 
tion length which is smaller than the pulse width. 

For reflections from the ocean, the most realistic case is a » o^ 

T f 

and » o, . The MSE for this case is 
T h 
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5 2 

+ S + ^ - 


(5.18) 


F^om this we see that the second-order terms are roughly proportional to 
— ^ . To prevent the second-order terms from becoming too large, should 
not be too small compared to a^. This means that the resolution of the 
receiver should not exceed the bandwidth of the returned signal. 

The ML estimator given in Eq. (5.3) correlates the received pulse with 
the logarithm of the mean pulse shape. It is interesting to consider 
correlating the received pulse at one wavelength with the logarithm of the 


received pulse at the other wavelength, i.e., 

^12 max-j^/ dt S 2 (t + x) in * 


(5.19) 


If Sj^ has a high signal-to-noise ratio, we have 


in Sj^(t) = in{S^(t) + ASj(t)} 

_ AS,(t) 

5 to S (t) +=^7-— 

^ S^(t) 


(5.20) 


Using Eq. (5.20), we find the bias and MSE of the estimator are given by 


Bias = 


/ dt S2(x+ S^(t) 

— 00 

oo 

/ dt S 2 (t + in S^(t) 


/ dt F(t) to F(t) 


/ dt F(t) to F(t) 


(5.21) 
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1 2 “* *“* 


MSE 


fa2 - 

— j j dt S2(t + t) to S.(t) 

> 

3t -00 ^ 

t=Ti2^ 


to S,(t.) + 


ASj(tj) 




to S, (t,) 




r^2' S^(t2) 


^l"”^2“'^12. 



fa2 - 

— - j dt S2(t + t) to S,(t) 

■*v 




3t -<*> 

L 

^^12j 


^ *^*^1 ^ ‘^‘^2 3t, 3t, ”^r^2 ■*■ ^2^1 


‘1 2 ''2 


T a T ® T 
1 2 12 


/ dt S^Ct + Tj2^ 


to S^(tj) to S^Ct^) + 


Cs^(tpt2> 


/ dt S2(t + Xj2) to Sj(t) 


(5.22) 


For a Gaussian return pulse shape, we have 




2 , °f 

®T <Nj> <N2> 


(<VS " <V'‘J (2.2 + ,2)3/2 


4a^<N^XN2> 4( o^ + oj)^ 


(5.23) 


where we have assumed to be large compared with and o^. Comparing 
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Eq. (5.23) with Eq. (5.16), we find correlating with the logarithm of the 

other pulse has slightly smaller first-order terras, but it has second-order 

4 4 

terms that are proportional to which can be much larger than those 

3 3 

of the correlation receiver, which are proportional to 

One can also find the peaks of both return pulses and use the separa- 
tion time between the peaks as the estimate, i.e.. 


T -2 = arg max (S 2 (t)} - arg max (S.(t)} 
t t 


(5.24) 


The estimator is unbiased because the mean .positions of the peaks are iden- 
tical at two wavelengths. The timing error can be found from Eq. (5.10) to 


MSE 


VaiJ^ Sj(t)i 


■2raax, 


t=t 


Imax. 


EJ^ S,(t)| 
df^ ^ 


t=t 


2maxJ 


EJ^ S,(t)| 
dt‘‘ ^ 


t=t 


ImaxJ 


3t^3t2 ^S2^^1*'^2^ 


t =t =t 5h^"2 

1 *^2max ^ ^ 


1 ^2 *^lmax , 


^2^*^2max^ 


^l^*^lraax^ 


(5.25) 

where t, and t_ are the locations of the peaks of S, and S„, respec- 
Imax 2max 1 2 

tively. For a Gaussian return pulse shape, this equals 


MSE 


(<N^> + <N2>) 


. 2 ^ 2 ^ 2,3 , 2 ^ 2, 
^°h~^ <^f + V (Of V 
3 2 2 2 3/2 
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2/ 2 ^ 2 ^ 2,3 
o_(a. + a, + o_) 


/ j_ . j_\ °h °f " 

U, kJ . 2 ^ 2.3/2 , 2 2 . , 2.3/2 

' 1 2 (oj^ + + 2c^) 


(5,26) 


For a target with negligible range spread so that o;^ = 0, the MSE reduces 


to 


- ( 4 ^ < 5 pr) 


2/ 2 ^ 2,3 

Of(a^ + op 
%, 1 . . 2,3/2 


(5.27) 


which is Independent of the speckle averaging. This is true because 
speckle does not distort the received waveforms and, therefore, does not 
change the positions of the peaks. Comparing Eq. (5.27) with Eq. (5.17), 
we find the peak-detection algorithm performs about the same as the corre- 
lation receiver for point targets. On the other hand, for a target with 


large-range spread, the performance of the peak detection algorithm is 


MSE 


(<Np <N2>) 2// o2 


2/2 ( 4 . 0 ^) 3 /^ 


(5.28) 


which when compared with Eq. (5.18) shows that peak detection is inferior 
to the correlation receiver. The difference is large when is large. 

This difference occurs because when is larger than and a^, shot noise 
and speckle will cause fluctuations within the received pulse, which make 
the peak position shift back and forth. 

We can also use the separation between the centroids of the two pulses 
as the estimate, i.e., 

Tj ^2 “ (Centroid of S 2 (t)) - (Centroid of S^t)) 


/ dt t S2(t) / dt t spt) 


/ dt S2(t) / dt spt) 


(5.29) 
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It can be shown that, under the assumption that and have the same 
waveform, the estimator is unbiased and has a MSE given by 


I J 


MSE 


•00 1 — • 00 — 00 / 

i + .■■ ■ . .... .■ — .1 ^ I . . ■ 

f \ * A 


<Nj>‘ 


<N^> 


/ 1 _1 . 1 _ . M 2 /_! _1 2 

V<Nj> <N2> Kj K27‘*r \<Nj> <N2>y®f 


(5.30) 


which is true for general waveforms if is taken to be the rms width of 
the time spread of the target. The MSE of the centroid algorithm depends 
only on the time spread of the received pulse while the performance of the 
correlation receiver depends on the slope of the received pulse. For a 
given time spread of Che received pulse, the correlation algorithm can be 
better or worse depending on Che bandwidth of the pulse. This will be 
further discussed in Section 5.4. 

In Table 5.1, we list the estimators considered above and their per- 
formances for the Gaussian mean pulse shape in the o_ » a, and ol » 

T f T 

case. Comparing the performances of these estimators for Che Gaussian 
'mean pulse shape, we find computing the centroids performs best. The 
correlation receiver follows. This is reasonable because for Gaussian 
pulse shape, Che implementation of the ML estimator shown in Eq. (5.3) 
reduces to calculating the centroids. This can be seen from the following 
derivations : 


ML 


arg max •< / dt S^(t) Zn S^(t + t^) 


arg max ■< / dt S^(t) Zn G( a, t + t^) 


L— <» 


TABLE 5.1 
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PERFORMANCES OF TIMING ALGORITHMS FOR GAUSSIAN MEAN PULSE SHAPE, ASSUMING 

>> a, AND >> a,. 

T h T f 


riming Algorithm 

Mean-Square Timing Error 



Correlation Algorithm 

.4.4 

Tjj “ arg max[/ dt 3j(t)S2<t+T) ] 



-1 -1 -1 -1 4 

Correlation Algorithm 
(pulse shape known a priori) 


- arg max [/ dt S 2 ( t)T 2 ( t+r) I 

1.54(<Sj>"^+k"^-KN2>"^+K2S <4 

- arg oax [/ dt S (c)T (t-*-t)] 

t -• ^ * 



(<Nj>'^+Kj^+<N2>"‘+K2^) 

Correlation with Log Waveform 

6 6 

^ m 

tj 2 • arg max [/ dt S 2 (t+T)la Sj(t) ] 

4o^<NjXS 2> 4(o^+op^jK, 

6 

+ (Kr*<N->"^+Kr‘<s,>"‘) — 

IZ »— 27? 

h f 

Correlation with log pulse shape (Optlaium) 
(pulse shape kaown a priori) 


• arg max [/ dt 32(t)taS2(t + t) ] 

(<S"‘>+K"^KS2>"‘+K2b 

- arg max [/ dt S^(t) InT^(fhT) ] 


Peak Detection 

- arg max lS 2 (t)l - arg max l3j(t)l 

5 

-1 -1 

(<N^> ‘ + <N,> ‘) i 

2/2o^ 

h 

5 


• (Centroid of Sj) ” (Centroid of sp 

m m 

/ dt t S (t) / dt t S (t) 

^ ^ ^1 -1 -1 -1 2 

- — ^ ^ «Nj> ‘ + Kj + <N,> ‘ + K, ) 0^ 

/ dt / dt 3^(t) 
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= arg max|/ dt[-(t + S^(t)/2o^]j> , 

which means x. satisfies 

/ dt(t + S.(t)/2a^|j ^ =“ 0 , 

i - |t^-x^ 

which is equivalent to calculating the centroid, 

09 

] dt tSj^(t) 


(5.31) 


(5.32) 


ML 


I dt S^(t) 


(5.33) 


However, for general waveforms, calculating the centroids could be far from 
optimal. 

5.3. Frequency Domain Representation of Timing Error 

From the above results, we restrict our attention to the correlation 
and centroid algorithms. For general waveforms, the performance of the 
correlation receiver has to be calculated numerically. The MSE can be 
expressed in terms of the bandwidth of the received signal. Equation (5.13) 
can be written as 


* - 3^ - 

/ dt S^(t) S2(t + x) 


T=X. 


J dm ^fs^(t)}^{S2(t + x^2>r 


12 


<N,XN„> » , , 

/ dm m |(|)^( oj) p 

— CO 


2tt 


(5.34) 


where 


-imt 


4) (m) = / dm F(t^ ® 

S 

= / d _pb-(_p,z) exp 


1 ^ 2 ,' 

— (Oj^ + Of) 


exp[imij^] 


(5.35) 
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'j)^(aj) is the Fourier transform of the mean received pulse shape. For any 
function a(t) we use^{a(t)} to denote its Fourier transform. In arriving 
at Eq. (5.34), we have made use of the following Fourier transform rela- 
tionship 

QD 00 

/ dt a(t)b(t) / d<i»^{a(t)}^{b(t)}* . (5.36) 

— 00 — oo 


Equation (5.12) can be written as 


OB 

I I <tj,t2)}^' 

4 TT —00 — oo 


r .2 
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where 




/ dt^ / dt 2 Rg (tj,t 2 ) e 1 1 e “ 2 


(5.38) 


R_ (o), , 01 .) is the two-dimensional Fourier transform of R (t, ,t„). In 
arriving at Eq. (5.37), we have made use of the following Fourier transform 
relationship, 

OO 00 00 oo ^ 

/ dt^ / dt2 a(tj,t2) b(t^,t2) = — 7 / . 4 ) }^{b(t^ ,t2) } 

— 00 —00 4 IT** — oo 

(5.39) 


which is the two-dimensional version of Eq. (5.36). Using Eqs. (5.34) and 
(5.37), we have 
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The Fourier transform can be calculated using the FFT algorithm and the 
integrations can be done digitally, although the computation time can be 
substantial due to the presence of two-dimensional transforms and integra- 
tions. However, by considering Eq. (5.40) alone, it is not clear what is 
the relationship between the timing error and the characteristics of the 
received signal. By substituting the explicit expressions for Rg, we can 
express the timing error in terms of the bandwidth of the received signal. 
The Intermediate steps are shown in Appendix B; the results are 
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B is the RMS bandwidth of the received pulses, and a is a dimensionless 
factor which typically varies between 1 and 2 , depending on the surface 
profile. For Gauss ian-pulse shape, a equals 1.54. Since the second-order 

3 

terras are roughly proportional to (l/a^^) , to keep the second-order terras 
from becoming very large, the resolution of the receiver (l/a^^) should not 
be too large. Equation (5.4) shows that timing error is inversely propor- 
tional to the signal bandwidth. This is expected since high-bandwidth 
signals will contain fine structures which improve the performance of the 
correlation timing algorithm. 

5.4. Comparison of the Correlation and Centroid Algorithms 

We will evaluate the performance of the correlation and centroid 
algorithms numerically for a set of simulated ocean return waveforms. The 
waveforms we use are a raised cosine superimposed with smaller scale modu- 
lations. That is 

F(t) = (1 + cos (27Tt/D)) + A( I - cos (2T7t/d))” , -D/2 ^ t ^ D/2 , (5.44) 

where d is smaller than D. Here the variable n is used to control the 
sharpness of the modulation. The larger the value of n, the sharper the 
peaks. This approximates the ocean return signal plotted in Fig. 3.5 for 
off-nadir pointing. 

We first compute the results for a smooth raised cosine, which is 
obtained by letting A equal zero. The rras timing error is plotted versus 
expected photocounts in Fig. 5.1 for a speckle number ~ ^2 * 50,000. 
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Figure 5.1. Comparison of the correlation and centroid algorithms for a 
simulated received pulse shape from a Gaussian ocean surface 
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The centroid algorithm performs better In this case, which Is expected 
because the centroid algorithm performs better for the Gaussian pulse shape 
and the raised cosine resembles the Gaussian. 

Next, we use a waveform that approximates the return from a sinusoidal 
ocean shown previously in Fig. 3.5. The approximation is done by letting n 
equal 4, A equal 0.94 and D/d equal 12. The results are shown in Fig. 5.2; 
the performance of the correlation receiver is significantly better, 
because the waveform has a relatively large bandwidth due to the presence 
of the sharp structures. The rras bandwidths of the waveforms shown in 
Figs. 5.1 and 5.2 are calculated to be 0.272 and 2.717 GHz, respectively. 

In Fig. 5.3, the power spectra of the two waveforms are plotted. Indeed, 
the waveforms shown in Fig. 5.2 have more high-frequency energy. While the 
performance of the centroid algorithm remains nearly the same for both 
waveforms, because the mean square widths of both waveforms are approxima- 
tely equal, the performance of the correlation algorithm improves signifi- 
cantly with the Increase of the signal bandwidth. From the above results, 
we conclude that, to achieve picosecond timing accuracy with ocean 
reflected pulses, the correlation algorithm should be used. In Chapter 6, 
we consider the effects of the footprint on the bandwidth of the received 
signal. 

In Fig. 5.4, we show the performance of the correlation receiver for 
three different values of speckle number. At the strong signal limit, the 
performance is limited by speckle. The amount of speckle averaging deter- 
mines the timing error. At the weak signal limit, the performance is 
limited by shot noise and the timing error is less sensitive to the number 
of speckle cells. In all cases, the timing accuracy improves as the signal 
strength and speckle number Increase. Single shot accuracies of ~10 psec 
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Figure 5.2. Comparison of Che correlation and centroid algorithms for a 

simulated received pulse shape from a sinusoidal ocean surface. 
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Figure 5. 3a. Power spectra of the mean received pulse shape shown In 
Figure 5,1. 
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Figure 5.3b 


Power spectra of the mean received pulse shape shown in 
Figure 5.2. 


RMS ERROR IN PS 


111 


PERFORMANCE OF CORREL.ATIOM RECEIVER 
1^. 100^. 10000. 10000EI 



100. 1000. 10000. 10000EI 


NUMBER OF PHOTOELECTRONS 


Figure 5.4. Effects of speckle averaging on the correlation algorithm. 
The mean pulse shape corresponds to Figure 5.2. 
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appear to be feasible with ocean surface targets when the expected photo- 
counts for both pulses approach a few thousand. The timing accuracy can be 
Improved further by averaging the measurements from many shots. 

5.5. Transmitter Effects 

In practice, the expected waveforms at the two wavelengths may not be 
exactly the same. This can be due to misalignment of the transmitted 
laser beams and different beam divergence due to the diffraction effects or 
due to turbulence. When the mean waveforms at the two wavelengths are not 
the same, the correlation receiver can become biased, with a bias given by 


Bias = - 
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/ dt T^Ct) T2(t + t^2) 
— 00 
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(5.45) 


and the MSE given previously by Eq. (5.40) becomes 

" * . . lour 

j dm^^ j dm^ m^m^ Rg^Cm^m^) Rs 2 (wpm 2 ) e ^ e ^ 

MSE = - — = . (5.46) 
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This result Is expanded Into a form suitable for numerical computation In 
Appendix C. As an example, suppose one channel has a very large footprint, 
so that It has a mean pulse shape which looks like Fig. 5.1. Suppose the 
other channel has a small footprint and a mean pulse shape that looks like 
Fig. 5.2. The bias of the correlation algorithm Is found to he zero 
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because although the waveforms are different they are not shifted with 
respect to each other. The timing error can be evaluated and the results 
are shown in Fig. 5.5. Compared with Figs. 5.1 and 5.2, we see a degrada- 
tion in the performance, which is expected, because the correlation between 
the two waveforms has decreased. 

5.6. Sampling Effects 

In practice, the output of the receiver or photodetector is usually 
sent into a waveform digitizer and the sampled and digitized data are 
recorded by a computer for future processing. If the waveforms are sampled 
at a fixed interval, t^, then the output of the waveform digitizer can be 
modeled as 

1 T 

o 

Y = / dt S(t) + n. ^ (5.47) 

(i-l)T ^ 

o 

and 

E{n^nj} = Var{n^} 6^^ , (5.48) 

where is the 1-th sampled data point, n^ is the additive noise 
introduced in the sampling and digitizing processes. Here n^ is due to the 
quantization in the digitizing process and the noise of the waveform digi- 
tizer. To avoid losses in the high frequency content of the signal, 
usually is chosen to be short compared to the Impulse response of the 

receiver. In this case, we have 

Y^ = S(iT^) . (5.49) 

Next, we consider the effects of n^ on various timing algorithms. 

Since the analyses of timing algorithms usually Involve differentiations, a 
discrete process like Y^^ can not be readily analyzed. To overcome this 
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Figure 5.5. Performance of the correlation algorithm when the mean pulse 
shapes at two wavelengths are different. One wavelength has 
mean pulse shape corresponding to that in Figure 5.1. The 
other wavelength has mean pulse shape corresponding to that 
in Figure 5.2. 
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difficulty, we caa approximate with a continuous process Y(t), 

Y(t) = S(t) + n(t) (5.50) 

and 


= Rg(tj,t2) + N^G(Oj.tj - t2) • ^5.51) 

where o^ is the correlation length of the noise, which is on the order of 
-the width of the .digitizer element spacing t ,N ,1s the noise energy. The 

ti + t2 ° ° 

function G(o 2 , g introduced to keep the total noise energy 

finite in the correlation process, and 02 is chosen to be large compared 
to the support of S(t), so that the noise power is essentially uniform. 

By substituting R^(tpt 2 ) with RY(tpt 2 > into the previous equations, 
we find additional error terms are Introduced. The MSB of the correlation 
receiver, given previously in Eq. (5.15), is increased by an additional 
amount MSE’ , which is given by 
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In Eq. (5.52), we have neglected the higher-order terras which are generally 
very small. By substituting the expression for into Eq. (5.52), we 
obtain 


MSB' 


1 


/2it M 


1 




(5.53) 


The computation leading to Eq. (53) is given in Appendix D. For a 
Gausslan-pulse shape, MSE' becomes 


MSE’ 
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<n^>7 


(5.54) 


From Eq. (5.54), we see that for the correlation receiver the additional 
timing error introduced by the sampling-digitizing process is usually 
small, because it is proportional to very large. 

For the peak detection algorithm, the additional error can be found 


from Eq. (5.25) to be 
3t^3t2 
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For a Gaussian received pulse shape, Eq. (5.55) becomes 
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which is proportional to N^/a^o^. large and the additional 
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error for the peak detection algorithm can be significant when the signal 

levels <Nj^> and <N^> are not very large. If peak detection is to be used, 

it is advisable to make the correlation length of the noise larger. 

can be done by smoothing or low-pass filtering the sampled sequence. 

For the centroid algorithm, we find the sampling-digitizing process 

Introduces an additional timing error which can be found from Eq. (5.30) to 

be given by 
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(5.57) 


This value is usually small for a reasonable signal strength. 

The effects of sampling on the timing accuracy for a Gaussian mean 
pulse shape are summarized in Table 5.2. We conclude that correlation and 
centroid algorithms are not sensitive to the noise Introduced by the 
sampling-digitizing process, but peak detection is more sensitive. For 
peak detection, it is advantageous to first smooth the sampled data. 

5.7. . ..Summary 

The ML estimator for differential arrival time cannot be Implemented 
due to the lack of prior knowledge of the mean pulse shape. Among the 
suboptimal estimators investigated, the correlation algorithm can achieve a 
signal shot accuracy of ~ 10 psec if the reflected pulse has time structure 
or glints. The performance of the correlation algorithm is found ti 
miprove with the bandwidth of the received pulse. 

If the mean received pulse shape at the two wavelengths is not the 
same, bias and additional error are Introduced into the correlation 
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TABLE 5.2. 

ADDITIONAL TIMING ERROR DUE TO SAMPLING NOISE 

Timing Algorithm Additional -Mean-Square Timing Error 

Correlation Algorithm (<N,>~^ + <N,>"^)2y^N ( of + oj + 

1 z o t n T c 

Centroid Algorithm (<N > ^ + <N-> N ( o^ - 0.25 a^) 

1 Z O Z 1 

«N^>"2 + <N2>’^) N^(a^ 


Peak Detection 
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algorithm. Therefore, the laser beams at the two wavelengths should be 
carefully aligned. 

The sampling of the receiver output can introduce additional noise. 

The correlation algorithm is found to be less sensitive to the sampling 
noise than the peak detection algorithm. 
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6. BANDWIDTH OF THE OCEAN REFLECTED PULSE 
From Chapter 5 we found that the timing accuracy of Che receiver 
depends on the bandwidth of the received signal. In this chapter, we study 
the effects of the surface profile and altimeter parameters on the band-, 
width of the received signal. The expected bandwidth of the received pulse 
that has been reflected by a Gaussian distributed ocean surface is eva- 
luated under the conditions that Che profile correlation length is much 
less or much greater chan the footprint radius. 



(6.3) 
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where R_(p, ,p„) is the normalized autocorrelation function of the surface 
c , —1 —2 

profile 


We first assume L » z tan 0^ and a quadratic autocorrelation function 


with given by 


» Po ) “ 1 ~ 


Lpi - ^2 


(6.4) 
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where L Is the correlation length of the surface profile. The case of 
L « z tan 0^ is treated later. 

;By . -substituting ,Eqs. (6.3) and (6.4) into (6.2) and making the 
following changes of variables 
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By substituting Eq. (6.7) Into Eq. (6.1), we obtain 
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For the case when L « z tan 9^, the autocorrelation function can be 


approximated by a delta function. Therefore, 
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In this case, the expected bandwidth Is 
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L « z tan 9_ 


Combining the results for both L » z tan 0^ and L << z tan 0^, we can 


write 
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(6.15) 


The results for a flat surface can be computed by letting a^ equal zero. 
When L » z tan 0^, If the received pulse shapes are dominated by beam cur- 
vature (tan 0^ is large) or surface effects ( is large), y will be small. 


If the curvature and surface effects are negligible, y will be large. 


On 


the other hand, when L « z tan 0^, y will be large when is large. In 
these two limiting cases, 4 simplifies to 


Y « 1 
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Y Y » 1 . 
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To evaluate the expression for the bandwidth and timing error, we 
2 2 1/2 

assumed 0 ( 0 ^^+ a^) « a^< z tan 0^. This means that the surface rms 

height is large compared to the point target response width but small com- 
pared to the radius of the footprint. For small footprints, the solution 
for L » z tan 9^ is probably the most applicable so that 


a((^ + o z tan 


MSB = 


■ <N2> Kj K 2 / 


cL 


L » z tan 




(6.17) 


For large footprints the solution for L « z tan 9^ is the most applicable, 
and 

, 2 

/I 1 1 1 \ 

’ L « z tan 9j, (6.18) 

where we have given only the first-order terms of Eq. (5.41) for simplicity. 

Because the wave height is typically of the order of tens of cen- 
timeters, a^c will usually be a few nanoseconds or larger. Picosecond 
time resolution is probably not feasible when the footprint is large. 
However, when the footprint is small and the surface correlation length is 
large, picosecond accuracy is possible if the system parameters are pro- 
perly chosen. The requirement for a large correlation length is equivalent 
to requiring a glint within the footprint. 
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7. LINK EQUATIONS 


In this chapter, we use the results obtained from previous chapters to 
estimate the expected timing accuracies of realistic system designs and to 
estimate the corresponding accuracies of the barometric measurements. 

The expected number of received signal photons was calculated in 
Ref. 4 for Gaussian ocean wave statistics: 


n „ „2 V 

- hT « L — 
o 


|r(o) 


z- 4it(S^ + 2 tan^ 0^) 


(7.1) 


n is the receiving system efficiency, hf^ is the signal photon energy, Q is 

2 

the total transmitted energy, and T^ is the two-way atmospheric transmit- 

2 

tance. The ocean reflectance [R(0)| depends on wavelength and water com- 
position. From the near IR to the near UV, the ocean reflectance is ~2 % 

2 

[22]. The factor (S 2 tan 8^) in Eq. (1) is related to the solid angle 

2 

of the reflected signal. In a perfectly calm sea S s 0, so the solid 

angle of the back-scattered signal is equal to the solid angle of the laser 

2 

beam. For an Isotropic surface, S is related to the MS surface roughness 
and profile correlation length (L): 


= 2 . (7.2) 

2 2 

Bufton et al. [23] recently reported measurements of = |R(0)| /4 ttS 
for a variety of sea states, nadir angles, and laser wavelengths. Their 
measurements were obtained at wavelengths of 337 nm, 532 nm, and 9.5 urn 
using an airborne lidar. At the visible and near UV wavelengths, the ocean 
reflectance (8^) for low wind speeds was ~6%. Equations (3.8) and (3.9) of 
Chapter 3 are consistent with these results. 

If we use this value for 8^ in Eq. (7.1) and the parameters listed in 
Table 7.1 for a satellite-based altimeter, the expected signal photocount 
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is <N> = 1300. The speckle number calculated using Eq. (3.17) is 
K = 50,000. In this case, speckle is negligible, and the shot noise is the 
dominant source of timing error. From the data plotted in Fig. 5.4, we see 
that the single shot timing error for this case is ~18 psec. The error can 
be reduced to ~2 psec by averaging the timing measurements from 100 shots. 
If the fundamental (1064 nm) and tripled (353 nm) YAG laser frequencies are 
used, a 2 psec timing error results in a pressure error of ~1 rabar (see 
Fig. 2.1). 


TABLE 7.1. 

PARAMETERS OF A SATELLITE-BASED LASER ALTIMETER 


Receiving System Efficiency : 10% 

2 

Receiver Aperture Area : 0.1 m 

Two-way Atmospheric Transmittance; 0.5 

Altitude of the Satellite : 400 km 

Laser Beam Divergence Angle : 100 prad 

Total Transmitted Energy : 250 mJ 
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8. AIRBORNE ALTIMETER EXPERIMENT DATA ANALYSIS 
8.1 . Introduction 

An aircraft experiment was conducted on September 7 and 9, 1983, at 
the NASA facility- on Wallops Island, VA. The experiment was to test the 
flight performance of the instrumentation and to measure the correlation 
properties of the raultlwavelength ocean reflected return pulses. GSFC per- 
sonnel including Abshlre, McGarry and Rowe conducted the flight experiment, 
while the data reduction and analysis were performed at the University of 
Illinois. 

The experiment was performed on board the NASA WFF Elect ra airplane. 

A picture of the airplane is shown in Fig. 8.1. The laser ranging system 
used in this experiment, including the laser, computer, waveform digitizer 
and other electronics, is the same equipment used in previous horizontal- 
path experiments [24], [25], [26]. Figure 8.2 is a picture of the computer 
system on board the aircraft. Figure 8.3 is a picture of the equipment 
rack on the aircraft. Figure 8.4 is a picture of the arrangement of the 
equipment on the aircraft. 

On September 7, before the flight, the altimeter was first tested by 
ranging horizontally to a cube corner reflector. The system worked well on 
the ground. However, after the aircraft took off and ascended to the 
ranging altitude, the receiver did not detect any reflected laser pulses. 
After the flight, the ranging system was carefully examined. The problem 
was caused by the vibrations that occurred during the flight which caused 
a misalignment between the receiving telescope and the detector. 

On September 3, before the aircraft took off, two ground calibration 
data sets data were collected by ranging to a white plate. Then the 
aircraft took off, headed east over the ocean, and the altimeter started 
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Figure 8.1. A picture of the NASA Electra airplane. 
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Figure 8.2. A picture of the computer system on board the airplane 
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Figure 8.3. A picture of the equipment rack on board the airplane 



131 



Figure 8.4. A picture of the interior of the airplane showing the laser 
power supply, the laser and the transmitting optics. 


132 

collecting data. During the three-hour flight, the system operated 
smoothly except for a temporary computer failure due to a loose circuit 
board. The problem was fixed during the flight. In this flight, data were 
taken over the ocean at both 305 m (1000 ft) and 1219 m (4000 ft) altitudes. 
They are the first data that furnish Information on ocean reflection pro- 
perties with 8 cm height resolution [27]. 

8.2. Data Analysis 

A typical recorded waveform collected at 305 m altitude Is shown In 
Fig. 8.5. The pulse on the left corresponds to the return at 355 nm (UV), 
and the pulse on the right Is Che return at 532 nm (Green). The time scale 
In the plot Is the number of digitizer elements. The temporal spacing 
between the elements Is 39 picoseconds so that Che full scale range Is 
20 nsec. The actual flight time of Che green pulse Is shorter, but It Is 
passed through an electrical delay line so that both pulses can be recorded 
on Che waveform digitizer side by side. Due to this arrangement, Che 
separation time between Che two pulses will decrease with altitude. Instead 
of Increase. 

On some of Che recorded waveforms, double peaks In both return pulses 
were observed. An example Is shown In Fig. 8.6. This occurred In some of 
the returns from Che ocean at both 305 m and 1219 m altitude and also In 
some of Che returns from Che white plate. The separation between the 
closely spaced two peaks exceeds 1 nsec In some returns. 

Because a separation of 1 nsec Is equivalent to a target range spread 
of 15 cm, which Is not Che case for the white plate target, and unlikely 
for Che ocean target. The laser has a beam divergence angle of approximately 
500 prad (full width at 10% points); therefore, the footprint has a size 
of ~15 cm at 305 m altitude and ~60 cm at 1219 m altitude. An ocean-wave 
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Figure 8.6. A waveform recorded at 305 m (1000 ft) altitude. 
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height difference of 15 cm within a 60 cm footprint is possible, but is 
much less likely within a 15 cm footprint. We conclude that double pulsing 
of the laser is probably the cause of the double peaks in the return pulse. 
In future experiments, the resolution of the receiver will be increased so 
that this ambiguity will no longer exist. 

After the flight, the collected data were first run through a sorting 
program, in which every waveform was scrutinized to reject unsatisfactory 
ones. A waveform is rejected if any of the following conditions occur: 
signal level too low (no energy detected at one or two wavelengths), signal 
level too large (outside the dynamic range of the waveform digitizer), and 
signal not recorded completely due to a triggering problem. The good wave- 
forms are then processed to extract the differential arrival time between 
the two pulses. 

The first technique used for timing was the correlation technique. 

This was done by first calculating the correlation coefficient of the two 

return pulses. The correlation coefficient p is defined as 

N 

I x(i) x(i + j) 

1=-1 


p(j) 


N -1 1/2 f j+N 

I I 


x(i) 


1/2 


( 8 . 1 ) 


U=1 J Li=j+1 J 

where x(l) represents the output of the I digitizer element and N is the 
number of points used in the correlation. p(j) is normalized so that if 
the two pulses have exactly the same shape it achieves its maximum value of 
one. We calculated the correlation coefficient instead of just the corre- 
lation function because the correlation coefficient also furnishes infor- 
mation on the degree of similarity of the pulse shapes at two wavelengths. 
The delay value that maximizes p was then used as the differential arrival 


time estimate. That is. 
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T,_ = [arg max p(j)] t , (8.2) 

12 j 

where is the temporal spacing of the digitizer elements. For this 
experiment is 39 picoseconds for all the data. Therefore, only 

assume quantized values in integer multiples of t^. In retrospect, the 
experiment could have been improved by using a faster digitizer sweep 
speed. However, at the time of the experiment, there was no way of knowing 
how wide the return pulse would be. Therefore, a slower sweep speed was 
used to make sure the entire return waveform could be recorded. 

The quantization problem can be partially remedied by fitting a curve 
through the correlation function and finding the peak of the fitted curve. 

A second-order polynomial was chosen as the fitting curve, because most 
return pulses resemble Gaussian pulses and a Gaussian is well approximated 
by a quadratic near its peak. We first tried using 5 or more points to 
perform the least square curve fitting. However, because the fitted curve 
did not pass through all data points, there were cases when the peak of the 
fitted curve was smaller than the original correlation peak. Therefore, we 
decided to use only 3 points around the correlation peaks for curve 
fitting. In this case, the fitted parabola will pass through all 3 points. 
The success of the 3-polnt curve fitting requires a high signal-to-noise 
ratio for the data points used. This requirement holds for the points near 
the correlation peaks. 

A history plot and histogram of the timing and the correlation coef- 
ficient are shown in Figs. 8.7 and 8.8 for the 305 m altitude data. 

In Fig. 8.9, two recorded waveforms collected at 1219 m altitude are 
shown. The timing results are shown in Fig. 8.10. Statistics of the 
correlation coefficient are shown in Fig. 8.11. 
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WAVEFORM NUMBER 

# OF POINTS GREATER THAN 8.448 ns = 0 

# OF POINTS SMALLER THAN 7.667 ns » 0 

TOTAL NUMBER OF WAVEFORMS = 255 

AVERAGE SEPARATION TIME ^ 204.6 ELEMENTS 

= 7.994 NANO SEC 

STANDARD DEVIATION = 1.612 ELEMENTS 

= 63 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.7. Timing history and timing histograms for the 305 m altitude 
data, using the correlation algorithm. Correlation coef- 
ficient is curve fitted around its peak. 




138 


h- 




u 

1.66 ■ 
6.99- 



o 




u. 

6.98- 

■ 


U- 




UJ 

6.97- 

• • . ■ 


o 


* 


o 

6.9€ 

u 

• 


z 

6.95' 

• 


o 




»- 

< 

6.94 



ID 

6.93^ 



cr 

CD 

6.92- 



o 

o 

6.91 

6.96- 

< -'■ 1 .■ 1 1 ' 1 1— .... t 

1 t 


> - I I 1- ( —— I 1 h——i \ 

1 36 72 i«7 143 178 213 248 284 326 H35 

WAVEFORM NUMBER 


# OF POINTS SMALLER THAN 0.900 = 5 


TOTAL NUMBER OF WAVEFORMS = 355 

MEAN » 0.985 STD = 2.3 E-02 


in 

K 

Z 

3 

O 

o 


cr 

UJ 

CD 

z> 


49 

36 

32 

28 

24 

26 

16 

12 

8 


4 t 


-M- 


, 1 1,1 > t 


iL 



6.96 6.91 6.92 6.93 6.94 6.95 6.96 6.97 6.98 6.99 l.( 

CORRELATION COEFFICIENT 


Figure 8.8. History plot and histogram of the correlation coefficient for 
the 305 m altitude data. Correlation coefficient is curve 
fitted around its peak. 
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Figure 8.10. Timing history and timing histogram for the 1219 m altitude 
data, using the correlation algorithm. Correlation coef- 
ficient is curve fitted around its peak. 
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WAVEFORM NUMBER 


# OF POINTS SMALLER THAN 0.900 = 0 


TOTAL NUMBER OF WAVEFORMS = 177 



CORRELATION COEFFICIENT 


Figure 8.11. History plot and histogram of the correlation coefficient for 
the 1219 m altitude data. Correlation coefficient is curve 
fitted around its peak. 
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For both 305 m and 1219 m altitude data, the correlation coefficient 
has average values approaching 0.99. This high average value of the 
correlation coefficient Indicates that the pulse shapes at the two wave- 
lengths are very similar even on the pulse-to-pulse basis. This is 
encouraging for the pressure measurement technique because similar wave- 
forms are essential. 

In Fig. 8.12, two recorded waveforms from a white plate target are 
shown. The corresponding timing and correlation coefficient results are 
given in Figs. 8.13 and 8.14. In Fig. 8.15, we plot two recorded waveforms 
from a cube corner reflector. The timing results are shown in Fig. 8.16 
and the results for the correlation coefficient are given in Fig. 8.17. 

In Table 8.1, we summarize the above results. From the results, we 
find the data obtained from ranging to the white plate have the smallest 
timing fluctuations, the highest mean correlation coefficient and also the 
smallest standard deviation of the correlation coefficient. It is reason- 
able that the returns from the white plate would have a higher signal-to- 
noise ratio than the returns from the ocean, but it seems that the returns 
from the cube corner reflector should have an even higher slgnal-to-noise 
ratio. This was not the case, probably because the ranging to the cube 
corner reflector was performed in daytime with the presence of a strong 
turbulence during the experiment. 

The data obtained at 305 m and 1219 ra altitudes can be used to verify 


the pressure measurement concept. From Appendix E, the differential arri- 


val time measured at 1219 m and 305 m should differ approximately by 


At = 1326(e 


-305/h -1219/h 


= 126.6 psec , 


(8.3) 


MAGNITUDE MAGNITUDE 


3 



TIME SCALE 



TIME SCALE 


. Two recorded waveforms obtained from ranging to a white 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.S2& ns = 0 

M OF POINTS SMALLER THAN 7.745 ns = 0 

TOTAL NUMBER OF WAVEFORMS = 349 

AVERAGE SEPARATION TIME = 207.1 ELEMENTS 

= 8.093 NANO SEC 

STANDARD DEVIATION » 1.030 aEMENTS 

» 40 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.13. Timing history and timing histogram for the data obtained 
from ranging to a white plate, using the correlation 
algorithm. Correlation coefficient is curve fitted around 
its peak. 
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# OF POINTS SMALLER THAN 0.900 » ‘ 0 


TOTAL NUMBER OF WAVEFORMS = 349 

MEAN =0.990 STD » 7.0 E-03 



CORRELATION COEFFICIENT 


Figure 8.14. History plot and histogram of the correlation coefficient 
for the data obtained from ranging to a white plate. 
Correlation coefficient is curve fitted around its peak. 
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Figure 8.16. Timing history and timing histogram for the data obtained 

from ranging to a cube corner reflector, using the correla- 
tion algorithm. Correlation coefficient is curve fitted 
around its peak. 
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TABLE 8.1. 

SUMMARY OF THE RESULTS OBTAINED FROM THE CORRELATION TECHNIQUE. THE 
90 % CONFIDENCE INTERVALS OF MEAN AND STANDARD DEVIATIONS ARE SHOWN. 


Data Set 

Separation 

Mean 

Time 

STD 

Correlation 

Mean 

‘Coefficient 

STD 

305m 

7995 psec 

63 psec 

0.986 

0.023 

altitude 





data 

[7989,8001] 

[60,67] 

[.984,. 988] 

(.022, .025] 

1219 m 

7888 psec 

91 psec 

0.989 

0.010 

altitude 





data 

(7877,78991 

[84,101] 

[.988, .990] 

(.009, .011] 

White 

8093 psec 

40 psec 

0.991 

0.007 

plate 





target 

[8089,8097] 

[38,43] 

(.990, .992] 

(.006,. 008] 

Cube 

8094 psec 

43.95 psec 

0.979 

0.008 

corner 





target 

[8090,8098] 

[41,47] 

(.978,. 980] 

(.007,. 009] 
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where is the atmospheric scale height, taken to be 8737 m. From the 
experimental data, the difference in differential arrival time is 106 
picoseconds. Taking into account the fluctuations in the aircraft altitude 
and the quantization in the waveform digitizer, the 20 psec difference bet- 
ween the theory and the experiment is considered to be within the accuracy 
of the experiment. 

For comparison, we processed the data collected at 305 m altitude and 
the data collected from ranging to a cube corner reflector using other 
estimation schemes. These two sets of data are chosen for comparison 
because the ocean is an example of a target with range spread, while the 
cube corner reflector is an example of a point target. 

From Chapter 5, we know one of the simplest techniques for estimating 
differential arrival time is to find the separation time between the peaks 
of the two pulses. For the 305 m altitude data, the results are shown in 
Fig. 8.18. The standard deviation of timing is 176 picoseconds, much 
larger than that of the correlation algorithm. For the cube corner reflec- 
tor data, the results are shown in Fig. 8.19. The peak detection yields a 
timing error of 49 picoseconds, which is not far from the 44 picoseconds of 
the correlation algorithm. This is consistent with the theory derived in 
Chapter 5, which states that, for range spread targets, the peak detection 
performs significantly worse than the correlation algorithm, while for 
point targets, the performances of the two estimators are close. Also, 
from Chapter 5, we know that the peak detection algorithm is more sensitive 
to noise introduced by the waveform digitizer than the correlation 
algorithm. This partially explains the slightly inferior performance of 
the .peak detection algorithm for the cube corner reflector data. 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.400 ns = 0 

# OF POINTS SMALLER THAN 7.600 ns = 4 


TOTAL WMBER OF WAVEFORMS = 355 

AVERAGE SEPARATION TIME = 203.A ELEMENTS 

» 7.956 NANO SEC 

STANDARD DEVIATION = 4.507 ELEMENTS 

= 176 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.18. Timing history and timing histogram for the data obtained at 
305 m altitude, using the separation between the two peaks 
as the differential arrival time estimate. 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.500 ns « 0 

# OF POINTS SMALLER THAN 7.700 ns * 0 

TOTAL NUMBER OF WAVEFORMS = 300 

A'^ERAGE SEPARATION TIME = 206.9 ELEMENTS 

- S.083 NANO SEC 

STANDARD DEVIATION = 1.257 ELEMENTS 

= 49 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.19. Timing history and timing histogram for the data obtained from 
ranging to a cube corner reflector, using the separation 
between the two peaks as the differential arrival time estimate. 
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The centroid detection algorithm was also used to calculate the dif- 
ferential delay. We used a total of 75 points around the peaks to calcu- 
late the centroids, and the timing results for the 305 ra altitude data are 
shown In Figure 8.20. Compared to peak detection, the timing error Is 
improved; It Is 118.7 picoseconds. The corresponding results for the cube 
corner reflector are shown in Fig. 8.21. The timing error Is 50 picosec- 
onds, which is slightly larger than the 49.1 picoseconds of the peak 
detection algorithm. This Is recognized to be caused by the inclusion of 
the ringing .noise of the .waveform digitizer In calculating the .centroids. 
The return pulse from the cube corner reflector occupied a width narrower 
than 75 points, so using 75 points In calculating the centroids Includes 
the undesirable data points that correspond to the ringing of the digi- 
tizer. We then used only 50 points in calculating the centroids. The 
results are shown In Fig. 8.22; the timing performance is improved to 
46 picoseconds. 

From Chapter 5, we know that the ML estimator correlates the received 
pulse with the logarithm of Its mean pulse shape. Because Che mean pulse 
shape Is not available, the tlL estimator can not be implemented. 

Correlating one return pulse with the logarithm of the other return pulse 
will not be optimum, but It Is Interesting to see how It performs. We 
first correlated the UV pulse with the logarithm of the green pulse. The 
results are shown In Fig. 8.23 for the 305 m altitude data. The standard 
deviation of timing is 84 picoseconds, which is not as good as the 63 pico- 
seconds of the correlation algorithm. In Fig. 8.24, two recorded waveforms 
and their logarithms are shown. Taking the logarithm enhances the small 
amplitude fluctuations, which are mostly due to ringing and digitizer 
noise. The quantization effects also become more apparent after taking the 
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Figure 8.20. Timing history and timing histogram for the data collected at 
305 m altitude. Results obtained by calculating the centroids 
of 75 points around the peaks. 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.526 ns = 0 

# OF POINTS SMALLER THAN 7.745 ns » 0 

TOTAL NUMBER OF WAVEFORMS » 300 

AVERAGE SEPARATION TIME = 206.8 ELEMENTS 

- 8.080 NANO SEC 

STANDARD DEVIATION * 1.291 ELEMENTS 

* 50 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.21. Timing history and timing histogram for the data collected 

from ranging to a cube corner reflector. Results obtained by 
calculating the centroids of 75 points around the peaks. 
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WAVEFORM NUMBER 

i OF POINTS GREATER THAN 8.S26 ns = 0 

# OF POINTS SMALLER THAN 7.745 ns = 0 

TOTAL ^JUMBER OF WAVEFORMS = 300 

A'^AGE SEPARATION TIME = 206.6 ELEMENTS 


s 8.073 NANO SEC 
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Figure 8.22. Timing history and timing histogram for the data collected 
from ranging to a cube corner reflector. Results obtained 
by calculating the centroids of 50 points around the peaks 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.370 ns = 0 

# OF POINTS SMALLER THAN 7.589 ns = 0 

TOTAL NUMBER OF WAVEFORMS = 355 

A'.'ERAGE SEPARATION TIME = 202.9 ELEMENTS 

= 7.928 NANO SEC 

S'ANDARD DEVIATION = 2.160 ELEMENTS 

= 84 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.23. Timing history and timing histogram for the data collected at 
305 ra altitude. Results obtained by correlating the UV pulse 
with the logarithm of the green pulse. 
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logarithm. The corresponding results for the cube corner reflector data 
are shown in Fig. 8.25; the timing error is a large 105 picoseconds. 

Looking at the timing histogram in Fig. 8.25, we find the poor timing per- 
formance is caused by the presence of a few estimates that are signifi- 
cantly below mean. Two recorded waveforms and their logarithms are shown 
in Fig. 8.26. Again, we find ringing becomes very pronounced after taking 
the logarithm. 

We then correlated the logarithm of the UV pulse with the green pulse; 
the results are -shown in Fig. 8.27 for the 305 m altitude data. The .timing 
error is 68 picoseconds, better than the previous case, but still not as 
good as the correlation algorithm. The results for the cube corner reflec- 
tor data are shown in Fig. 8.28; the timing error is reduced to only 
43 picoseconds, which is close to the 44 picoseconds of the correlation 
algorithm. This can be explained by noting the fact that, for a point 
target, the individual received pulse has the same shape as the mean pulse 
shape, which is just the impulse response of the system. The difference is 
that the Individual return is corrupted by ringing and noise of the digi- 
tizer. Therefore, for -a point -target , if -ringing is. not severe, corre- 
lating one received pulse with the logarithm of the other pulse is a good 
approximation of the ML estimator. This explains why it performs quite 
well. For both data sets, we find the algorithm that takes the logarithm 
of the UV pulse performs better than those taking the log of the green 
pulse, probably because the UV pulse has less of a ringing problem than the 
green pulse. 

Finally, we correlated the logarithms of both received pulses. The 
results for the 305 m altitude data are shown in Fig. 29; the timing error 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.448 ns = 0 

« OF POINTS SMALLER THAN 7.467 ns = 0 

TOTAL NUMBER OF WAVEFORMS = 300 

AVERAGE SEPARATION TIME = 204.9 ELEMENTS 

* 8.007 NANO SEC 

STANDARD DEVIATION = 2.698 ELEMENTS 

« 105 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.25. Timing history and timing histogram for the data collected 
from ranging to a cube corner reflector. Results obtained 
by correlating the UV pulse with the logarithm of the green 
pulse. 
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WAVEFORM NUMBER 


# OF POINTS GREATER T>«N 8.408999 ns = 0 

« OF POINTS SMALLER THAN 7.628000 ns = 0 

TOTAL NUMBER OF WAVEFORMS = 355 

A'.ERAGE SEPARATION TIME = 204.0896 ELEMENTS 

= 7.972250 NANO SEC 

STANDARD DEVIATION = 1.764449 ELEMENTS 

= 68.92381 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.27. Timing history and timing histogram for the data collected at 
305 m altitude. Results obtained by correlating the logarithm 
of the UV pulse with the green pulse. 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.526 ns = 0 

» OF POINTS SMALLER THAN 7.745 ns = 0 

TOTAL NUMBER OF WAVEFORMS = 300 

AVERAGE SEPARATION TIME * 206.9 ELEMENTS 

= 8.084 NANO SEC 

STANDARD DEVIATION = 1.098 ELEMENTS 

= 42 PICO SEC 



SEPARATION TIME (rrs) 


Figure 8.28. Timing history and timing histogram for the data collected 
from ranging to a cube corner reflector. Results obtained 
by correlating the logarithm of the IJV pulse with the green 
pulse. 
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WAVEFORM NUMBER 


# OF POINTS GREATER THAN 8.331 ns » 0 

# OF POINTS SMALLER THAN 7.550 ns = 0 

TOTAL. NUMBER OF WAVEFORMS = 355 

AVERAGE SEPARATION TIME = 201.6 ELEMENTS 

= 7.879 NANO SEC 

STANDARD DEVIATION = 2.651 ELEMENTS 

= 103 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.29. Timing history and timing histogram for the data collected 
at 305 m altitude. Results obtained by correlating the 
logarithms of both pulses. 
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Is 104 picoseconds. The corresponding results for the cube corner reflec- 
tor are shown In Fig. 8.30; the timing error Is 100 picoseconds. After 
taking the logarithms, the correlation between the two pulses decreases, 
because the uncorrelated noises are emphasized by the logarithm. From the 
above results, we conclude that the logarithm algorithms should be applied 
only to mean waveforms, or signals with very high slgnal-to-nolse ratios. 

The results from the different algorithms are summarized In Table 8.2. 

8.3. Altitude Dependence of the Received Signal Energy 

During the experiment. It -was found that the return signal level 

remained about the same when the altitude of the aircraft was Increased 

from 305 m (1000 ft) to 1219 m (4000 ft). According to the link equation, 

-2 

the signal level should vary Inversely with the altitude squared (z ). 
Therefore, the signal level at the 305 m altitude should be 16 times the 
signal level at 1219 m. 

-2 

This deviation of the signal level in Ildar systems from z dependence 
has been studied by Harms [281, [29]. He found there are two mechanisms 
that cause power loss at sraall-to-medium target distances. First, in a 
system with shadowing of the main .mirror by a central obstruction, such as 
a secondary mirror, the overlap of the transmitted beam with the receiver 
field of view is often Incomplete. Second, because the telescope is 
usually focused at Infinity, the light backscattered from a target at 
sraall-to-medlum distances Is not completely focused onto the detector. 

Both .effects cause power losses that can be significant at small target 
distances. 

In Appendix F, the received signal level is derived taking into 
account both effects. In Fig. 8.31, the relative signal level is plotted 
versus altitude for four different detector sizes. The detector Is assumed 
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WAVEFORM NUMBER 


# OF POINTS greater THAN 9.408 ns = 0 

# or POINTS SMALLER THAN 7.628 ns - 0 

TOTAL NUMBER OF WAVEFORMS » 300 

AVIRAGE SEPARATION TIME = 203.8 ELEMENTS 

* 7.963 NANO SEC 

STANDARD DEVIATION » 2.563 ELEMENTS 

« 100 PICO SEC 



SEPARATION TIME (ns) 


Figure 8.30. Timing history and timing histogram for the data collected 

from ranging to a cube corner reflector. Results obtained by 
correlating the logarithms of both pulses. 
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TABLE 8.2. 

COMPARISON OF THE RESULTS FROM DIFFERENT TIMING ALGORITHMS. THE 90 % 
CONFIDENCE INTERVALS OF THE TIMING ERRORS OF THE ALGORITHMS ARE SHOWN. 


Algorithm 


Timing Accuracy 


305 m Altitude Data 


Cube Corner Data 


Correlation 

63 psec 

44 psec 


[60,67] 

[41,47] 

Peak detection 

176 psec 

49 psec 


[166,188] 

[46,53] 

Calculate 

118 psec 

46 psec 

centroid 

[111,126] 

[44,50] 

Correlate the 

UV pulse with 

84 psec 

105 psec 

log of the 
green pulse 

[80,90] 

[99,114] 

Correlate log 

of the UV pulse 

69 psec 

43 psec 

with the green 
pulse 

[65,74] 

[40,46] 

Correlate log 

104 psec 

100 psec 

of both pulses 

[98,111] 

[94,108] 
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Figure 


RMS Laser Beam Width =0.117 mrad 
Telescope Focal Length = 14.47m (570 inch) 
Primary Lens Radius = 20.32 cm (8 inch) 
Central Obstruction Rodius =5.08 cm (2 inch) 
Telescope is Focused at Infinity 



8.31. Received signal level versus altitude for a telescope with 
central obstruction. 
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to be at the focal plane, or the telescope is focused at infinity. From 

-2 

the figure we see that, if the detector is very large, the z dependence 
is maintained, because all the energy collected by the telescope is 
detected and the assumption behind the link equation is valid. For the 
aircraft experiment, the detector is situated at the focal point of a 
telescope with an effective focal length of 14.47 m. The detector of 
Interest has a 3.175 mm (1/8 inch) radius. From Fig. 8.31, we see the 
theory does predict a nearly equal signal strength at 305 m (1000 ft) and 
1219 m (4000 ft) altitudes. When calculating the -maximum allowable alti- 
tude based on ground-based tests over short paths, the effects must be 
taken into account. 

In Fig. 8.32 we show the results for the case of no central obstruc- 
tion. In this case, light loss at close distances is due to the fact that 
the target is not focused onto the detector. The out-of-focus Intensity 
pattern spreads out over a larger area, and only a portion of it is sensed 
by the detector. In Fig. 8.33, we show the effects of the telescope 
focusing on the relative signal level. In the figures, when the telescope 
is focused at 305 m, the detector is the position where the image of a 
target 305 m away will form. In the figure, we see that the maximum 
received energy is obtained at a given altitude when Che target Is focused 
onto the detector. As a comparison, in Fig. 8.34, we plot the results for 
a telescope without the central obstruction. 

From the above results, we see that the altitude dependence of the 
received signal energy observed in the experiment can be accounted for by 
the central obstruction of the telescope and Che finite size of the detec- 
tor. For future experiments, these results can be used to more accurately 
predict the received signal level. 
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0 
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32. Received signal level versus altitude for a telescope without 
central obstruction. 
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RMS Laser Beam Radius = 0.1 17 mrad 
Telescape Facal Length = 14.47m (570 inch) 
Primary Lens Radius = 20.32 cm (8 inch) 
Central Obstructian Radius =5.08 cm (2 inch) 
Detectcr Radius =3.175 mm (0.125 inch) 



ALTITUDE IN METERS 


Figure 8.33. Effects of telescope focusing on the received signal level. 
Telescope has central obstruction. 
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RMS Laser Beam Radius = 0.117 mrad 
Telescope Focal Length = 14.47m (570 inch) 
Primary Lens Radius = 20. 32cm (8 inch) 
Detector Radius = 3.175 mm(O.I25 inch) 



Figure 8.34. Effects of telescope focusing on the received signal level. 
Telescope has no central obstruction. 
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3.4, Summary 

In this chapter, we analyzed the data obtained from the first airborne 
multiwavelength altimeter experiment. The data were processed with dif- 
ferent timing algorithms. Taking into account the confidence intervals of 
the results obtained, we found that the correlation algorithm probably has 
the best timing accuracy. 

A study was done on the altitude dependence of the received signal 
level. The results obtained explained the signal level observed during the 
experiment. 
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9. CONCLUSIONS 

Atmospheric pressure can be measured from a satellite by observing the 
change with frequency In the optical path length from the satellite to the 
earth's surface. Because the optical path-length difference Is relatively 
Insensitive to surface temperature and water vapor pressure, only rough 
estimates of these parameters are required. Therefore, the major error 
source Is likely to be the differential path-length measurement. For a 
nadir viewing altimeter using the fundamental and tripled YAG laser fre- 
quencies, the differential propagation delay Is approximately 52 cm. For a 
pressure accuracy of 1 mbar the differential delay must be measured with an 
accuracy of 0.3 mm or 2 psec. Instrumentation such as the streak tube 
camera which Is capable of meeting these timing requirements Is currently 
available. 

To achieve the high timing accuracy needed, the receiver has to cope 
with not only shot noise but also time-resolved speckle. Time-resolved 
speckle causes fluctuation of the received energy, distorts the received 
pulse shape and degrades the timing performance. This study showed that 
generally the mean square timing error Is Increased by a factor of <N>/K, 
due to the presence of time-resolved speckle. Fortunately, for realistic 
altimeter systems, K is usually very large (~50,000). 

Among the timing algorithms considered, we found the correlation 
algorithm to be most promising. The accuracy of the correlation technique 
improves as the bandwidth of the received signal increases. The success of 
the technique depends on the presence of specular reflections or "glints" 
within the footprint which preserve the high-frequency content of the 
transmitted pulses. Under this condition, we found that, by using Che 
correlation algorithm to determine the differential propagation delay, 
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single shot accuracies of ~10 picoseconds can be achieved when the expected 
signal photocount approaches a few thousand. 

As for the surface profile, we found that for glints to occur within 
the footprint, the profile should have a long correlation length, so there 
will be large areas where the surface elevation reioalns nearly constant. 

The condition is satisfied by sinusoidal and trochoidal waves. 

Experimental results from the first airborne altimeter experiment con- 
ducted at Goddard's Wallop Flight Facility show a high degree of correla- 
tion between the received pulse shapes at the visible and UV wavelengths. 
This is essential to the success of the pressure measurement technique. 
Although this initial experiment was not designed to provide data of suf- 
ficient accuracy to pemlit actual barometric measurements, they were used 
to verify the concept of the pressure measurement technique. We processed 
the data with various timing algorithms. The results show a general 
agreement with the theoretical performances of these algorithms. The 
results also indicate the superiority of the correlation technique. 

Currently, work is underway at NASA-Goddard to upgrade the altimeter 
system. An image intensified streak tube camera with timing resolution of 
~2 picoseconds will be used to replace the waveform digitizer used in the 
initial experiments. Extensive airborne experiments using the new system 
are scheduled for the summer of 1984. With the much higher timing accuracy 
of the new system, the data collected will be used to demonstrate the 
feasibility of the technique for the pressure measurement over the ocean. 

This dissertation work began with the study of the statistics and 
waveforms of the ocean reflected laser pulses. The results are not only 
essential to the pressure measurement technique, but also have applications 


176 


in remote sensing of sea states. Then, the statistical properties of time- 
resolved speckle and estimation of arrival times in the presence of both 
shot noise and speckle vere studied. The study was done in the general 
context of laser altimetry and the results obtained are new in the area. 

For completeness, an example of partially developed time-resolved speckle 
was also studied. 

The central problem of the dissertation work is the estimation of dif- 
ferential arrival time. Various algorithms were proposed and studied. The 
relation between timing accuracy and bandwidth of the received signal was 
established. These results enable the prediction of the accuracy of a 
realistic pressure measurement system. 

Data analysis was also a major part of the dissertation work. The 
data were processed with different timing algorithms and the results were 


compared with the theory 
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APPENDIX A 

EVALUATION OF THE PERFORl-IANCE OF THE COHRELATION ALGORITHM 
FOR A GAUSSIAIT MEAN RECEIVED PULSE SHAPE 

From Chapter 4, the mean and covariance functions of the received 

signal can be written explicitly as 


Sj(t) = <Nj> G t) . 


(A.l) 


S2(t) = <N2> g(/oJ + <^ + q2, t - Tj2> 


(A.2) 


Cs = <Nj> G(/2a^, t^ - t2> G(/o^/2 + + o^. (t^ + t^)/2) 


+ — G(/2/oJ + o^, t^ - t2> G(/<^/2 + 0^/2 + (t^ + t2)/2) 
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(A. 4) 


The denominator on the right-hand side of Eq. (5.15) can be calculated 
by making use of Eqs. (A.l) - (A. 4), 
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2 2 
<N^>'‘<N2> 


16ir(o^ + oj + 0^)^ 


The numerators of the four terms on the right-hand side of Eq. 
are evaluated as follows: 
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Following the same derivations, we have 
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Finally, by using Eqs. (A.5 - A.9), we obtain 
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The results are used in Chapter 5 
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APPENDIX B 

DERIVATIONS LEADING TO THE FREQUENCY DOMAIN REPRESENTATION OF TIMING ERROR 


From Chapter 4, for a general mean waveform F(t), we have 


S^(t) =■ <N^> F(t) , 
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where we have assumed the width of F(t) to be much larger than and a^. 
The mean square error Is written here for convenience: 
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The numerator on the right-hand side is calculated as follows; 
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By the same reasoning, we have 
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The other non-zero term is 
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By substituting Eqs. (B.6-B.9) into Eq . (B.5), we obtain Eq. 


(B.9) 

(5.41). 
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APPENDIX C 


EXPANSION OF EQUATION (5.46) 


Following the same derivations shown In Appendix B, and making changes 
for the difference In mean waveforms, we have 
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This result Is used In Chapter 5. 
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APPENDIX D 

DERIVATIONS LEADING TO EQUATION (5.53) 


>2 « 


^ ^ lOB — O 


/ dt S, (t) 

-» ^ 3 t^ 


S 2 (t + t) 


T=T 


12 


SjCti + tj) Sj(tj * Tj) R„(t,,tJ 


’^r'^2^’^12 


/ dt s,(t) -^• 
-» ‘ 3t 




T=T 


12 


C, + t 


N^<<Ni> + <N^> ) / dtj / dt2 F(tj) F(t2> G(a^, tj - t2> G 02 * 


1 2 


I ■* ~ "I 2 

<N^>^<N2>^<{ / dt F(t) F(t)| 


N^(<Nj>^ + <N2>^) / dt F^(t) G(02,t) 


<N^>^<N2>^|/ dt F(t) F(t)| 

C9 

No«Ni>^ + <N2>^)G(02,0) I7 / u^^gCw)!' 


12 


<N, 


1>^^V 117 / do. (D^ 1 *^( 0 ,) I 


N /ziT 

O 




1 


/ do) ( 0 ^ |<^ (la) 1^ 

*49 


(D.l) 


where we have made use of the fact that Oj is very small and 02 is large 


compared to the width of the received waveform 
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APPENDIX E 

EXPECTED DIFFERENTIAL DELAY FOR THE AIRCRAFT EXPERIMENT 


For Che aircraft experiment, we can write Che atmopsherlc correction 


as 


*^alrcrafC 10 ^ N 


AC 2 2 / 
■'r 


sin E 


^dr . 


(E.l) 


ocean 

where we have neglected the term due to the difference between Che 
geometric lengths of the ray and straight-line path. 

From Ref. [30], the group refractlvlty Is given by 


N 2 80.343 f(X) “ . 

s * 


(E.2) 


If a zero temperature lapse rate Is assumed [31] the pressure P will vary 
exponentially with altitude: 


where 


Pg exp 


RT 


h 


s J 
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s MG 


(E.4) 


h^ Is the atmospheric scale height and 


Pg = surface pressure, 

h =■ altitude of the aircraft, 

M = molecular weight of air, 

R =» universal gas constant, 

T = surface temperature (®K), 
s 

G = acceleration of gravity. 
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By substituting Eqs (E.2) and (E.3) into Eq. (E.4), we have 
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where 


P 

N = 80.343 f(X) . (E.6) 

S T 

S 

For E “ 90®, P = 1010 mb and T = 300*K, the scale height equals 8787 m 
and the difference of the atmopsheric correction at two wavelengths is 


AAC 5 4.76(m)(l - exp(-h/hg] )(f(Xj) - f(X 2 >) . (E.7) 

For the aircraft experiment, X^ = 0.532 pm and X 2 » 0.355 pm which results 
in 


AAC “ 0,3978(m) (1 - exp(-h/h ]) 

s 

= 1326 (ps) (1 - exp{-h/h^]) 


(E.8) 


The results are used in Chapter 8 
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APPENDIX F 

EFFECTS OF TELESCOPE CENTRAL OBSTRUCTIONS ON RECEIVED SIGNAL LEVEL 

To study the effects of the central obstruction of the telescope on 
the received signal level, we assume a target with uniform reflectivity and 
negligible range spread. The analyses here closely follow that of Refs. 
[28] and [29] . 

Let tp and r^^ be the radius of the primary mirror and central obstruc- 
tion (secondary mirror), respectively. Referring to Fig. F.l, the detector 
is situated on the optical axis, at a distance D behind the primary lens, 
while the obstruction or secondary mirror of the telescope is located at a 
distance S^ in front of the lens. The primary lens forms an image of the 
target on to the image plane. For a target at distance z, if the lens has 
a focal length f, the image is located at distance I(z) behind the lens, 
where I(z) is related to f and z by the thin lens formula: 

1/I(z) + 1/z = 1/f . (F.l) 

The light back scattered from a point r^ of the target is focused onto 

the image plane and causes an irradlance dS(r^,z) at point r^, where 
l( 2 ) 

r. => r . A pencil of light back scattered to a surface element dA 

X 0 2 

around r^^ will illuminate an annulus in the detector plane around r^ with 
an irradlance of 


dS( r^,z)dA 



dS(r^,z) 


(Fi2) 
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Figure F.l. Geon®try of the optical setup. 
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where S(r^,z) and SCr^^.z) are the Irradlances in the detector and image 

plane, respectively, r^ is related to r^ by r^ = Y ^zT ^i* ^c outer 

radius of the annulus, which is related to r by 

P 


= r - P 

’^c p I(z) 


(F.3) 


Equation (F.3) can be obtained from observing the geometry shown in Fig. 
F.2. Equation (F.3) is still valid when D > I(z) if we repalce (I(z) - D) 
by its absolute value. 

The total irradiance at r^ is provided by all the pencils of rays for 
which the principal rays Intersect the detector plane within a circle 
around r^ with radius r^ excluding those rays within a circle with radius 
r — . The corresponding marginal points in the image plane form an annu- 
lus with an outer radius r^ given by 


I(z) 

r = r — r — , 

m c D 


(F.4) 


and an inner radius equal to r — This can be seen from the geometry 

P 

shown in Fig. F.3. 

Assuming a Gaussian footprint with a radius at exp(-l/2) point equal 
to z tan 9^, the irradiance at the image plane is given by 
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2Tra^ 


(F.5) 


where 


a 


I 


z tan 9_ 
T 


I(z) 

Z 


I(z) tan 0^ . 


(F.6) 


Since we are only interested in the relative signal level, in Eq. (F.5) 
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Figure F.2. Ray tracing of the telescope. 
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we have omitted a proportional constant for simplicity. The constant takes 
Into account transmitted laser power, reflectivity of the target and 
atmospheric transmittance. 

By making use of Eqs. (F.5) and (F.2), we can express the Irradlance 
at the detector plane as 


S(rj^.z) 


2/. A(r^.|r^)/A(r_^-j^. ir^) 


irr_ 1 - 




2o: 


— 2 

z 


2iro‘ 


(F.7) 


X ** b I 

where A(r^,-|g- r^ )/A(r^ denotes the Integration region, an annulus 

P I 

with outer radius r , and inner radius r — centered around -;r r 

ra m r D d 

P 

By making changes of variables, we can write 


S(r^,z) 


1 


I 


r^z^ 2iro^ \(r^,0)/A(r^ O) 

_2 2 P 

I r. 


2 2 2 


2 2 

2o,D r 

e ^ / ” r^ dr r e 

•^m — 


d r exp 


2 a: 


I I |2‘ 

pJjl 
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(F.8) 


where I is the modified Bessel function of the first kind, zeroth order, 
o 

The total power received by the detector E(z) can be obtained by 
Integrating S(r^,z) over the detector area. 
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E(z) = 2 tt / dr^ SCr^.z) ^ ^ 


2tc R 


2ah^ r 


E_ f d 

2 2 2-' “‘d ‘d 


TqZ Oj o 


m r 


2 0j /Irr^^ 
re I r- 

°Ky 


Again by making changes of variables, we can write E(z) as 

2 

r 


E(z) 


2ff 




4 2 „ 

z tan 0„ 


z 

R:. TT 
f d D . 

J dr r e 
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2 2 
2z tan 0„ 


z z 

“ 7 ■*■ d"! 


2 2 
2z^ tan^ 0„ 


dp p e 


pr 
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\z tan 0_ 


(F.9) 


(F.IO) 


If the radius of the detector Is very large, l.e., R^ » z tan 0^, 
Eq. (F.IO) can be simplified to 


E(z) 


2v 




-f-f) 


2 2 
Z 


dp p 


<r2 - 2) 

^ p b-* 


(F.ll) 


which Is Inversely prportlonal to the distance squared. This functional 
dependence Is expected when all the energy collected by the telescope is 
Intercepted by the detector. On the other hand, if the detector Is very 
small so that R^ « z tan 0^, In Eq. (F.IO) can be approximated by 1 and 


E(z) simplifies to 




(F.12) 


When Che deCecCor Is situated at Che image plane, i.e., when the 


target is focused onto the detector, Che power received by the detector is 



(P.13) 


which is generally larger than Che power obtained when the detector Is at 
any other position. 
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